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A new model realisation of the vector curvaton paradigm is presented and analysed. The model 
consists of a single massive Abelian vector field, with a Maxwell type kinetic term. By assuming 
that the kinetic function and the mass of the vector field are appropriately varying during inflation, 
it is shown that a scale invariant spectrum of superhorizon perturbations can be generated. These 
perturbations can contribute to the curvature perturbation of the Universe. If the vector field 
remains light at the end of inflation it is found that it can generate substantial statistical anisotropy 
in the spectrum and bispectrum of the curvature perturbation. In this case the non-Gaussianity 
in the curvature perturbation is predominantly anisotropic, which will be a testable prediction in 
the near future. If, on the other hand, the vector field is heavy at the end of inflation then it is 
demonstrated that particle production is approximately isotropic and the vector field alone can give 
rise to the curvature perturbation, without directly involving any fundamental scalar field. The 
parameter space for both possibilities is shown to be substantial. Finally, toy-models are presented 
which show that the desired variation of the mass and kinetic function of the vector field can be 
realistically obtained, without unnatural tunings, in the context of supergravity or superstrings. 
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I. INTRODUCTION 



Cosmic inflation is the most compelling solution to the horizon and flatness problems of the standard, hot big 
bang cosmology. As an additional bonus, inflation can also generate the primordial curvature perturbation £ of 
the Universe, which is necessary for the formation of the observed structures such as galaxies and galactic clusters. 
Fortunately, a snapshot of the primordial curvature perturbation is reflected onto the Cosmic Microwave Background 
. (CMB) radiation, whose temperature perturbations can provide detailed information about (, such as the amplitude 
' and scale dependence of its spectrum. The importance of this information lies in its discriminatory power as it offers, 
qq , so far, the only tool which allows the testing of individual inflation models. Detailed cosmological observations, not 
H ' only of the CMB but also of the distribution of structures in the Universe (e.g. galaxy surveys), are in agreement 

j. | | with the generic predictions of inflation for Indeed, the observational evidence supports the existence of a nearly 

scale-invariant spectrum of predominantly adiabatic and Gaussian superhorizon perturbations. Inflation typically 
' produces such perturbations through the amplification of the quantum fluctuations of suitable fields during a period 
of accelerated expansion of space. The process is called particle production and it is a manifestation of the Hawking 
' radiation process in de Sitter space. Now that observations of the products of this process are becoming ever more 
. , precise, model-builders are forced to construct inflationary models which are better and more closely connected with 
■ the theory, more detailed and rigorous, and more testable and falsifiable by comparison with the data. The question 
H | to be asked then is, which are those "suitable fields", whose quantum fluctuations are the ultimate sources of the 
- - i curvature perturbation and all the observed structures in the Universe? 

A very useful tool in theoretical physics and, by extension, in cosmology has been the fundamental scalar quantum 
field. Such scalar fields are ubiquitous in theories beyond the standard model such as supersymmetry (the scalar 
partners of the observed fcrmion fields) or string theory (string moduli fields). In cosmology, a scalar field is usually 
assumed to drive the dynamics of inflation (inflaton field) or even give rise to the accelerated expansion at present 
(quintessence field). It is not surprising, therefore, that the particle production process, responsible for the generation 
of C, has until recently been considered only in the context of scalar fields. So useful a tool is a scalar field that one 
tends to forget that, as yet, it is no more than a theoretical construction, whose physical existence is only conjectured, 
since no fundamental scalar quantum field has ever been observed. One could wonder therefore whether other kinds 
of fields contributed to the generation of the curvature perturbation in the Universe. What if scalar fields did not 
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exist? Could the galaxies form without them? 

Recently, the contribution of vector boson fields to the curvature perturbation is attracting growing attention. In 
the pioneering work in Ref. pj a mechanism was introduced, which could potentially allow a single vector field to 
generate the observed curvature perturbation. The mechanism was the first such suggestion in the literature and was 
based on the curvaton idea, which was developed for scalar fields in Ref. Q (see also Ref. Q) in order to alleviate 
the fine-tunning problems of inflation model-building The vector curvaton mechanism assumes that the vector 
field has a negligible contribution to the energy budget of the Universe during inflation, when on the one hand it 
becomes homogenised, but on the other hand can also obtain a scale-invariant superhorizon spectrum of perturbations 
due to some appropriate breaking of its conformality. After inflation the density parameter of the vector field grows, 
especially after the mass of the vector field becomes important (larger than the Hubble scale). When this happens, 
the zero-mode undergoes coherent oscillations during which the vector field acts as pressureless, isotropic matter 
Hence, it can dominate (or nearly dominate) the Universe without generating excessive anisotropic expansion, 
avoiding thus the primary obstacle for using vector fields to generate C- When the field dominates the Universe it 
imposes its curvature perturbation spectrum according to the curvaton scenario • 

The vector curvaton mechanism is a paradigm in search of a model, much like inflation itself. The main ingredient 
necessary is a mechanism for the generation of a flat spectrum of vector field perturbations during inflation. This is 
a problem which was originally addressed by efforts to generate a coherent primordial magnetic field during inflation, 
in order to source the magnetic fields of the galaxies (for a review see Ref. Q). As realised in Ref. [f|, the particle 
production process can operate with vector fields only when their conformal invariance is appropriately broken. A 
number of proposals have been put forward for this to occur during inflation ■ For example, as discussed in 
Ref. a flat spectrum of vector field perturbations is generated if the effective mass-squared of the vector field 
during inflation is m 2 s = —2H 2 , where is the inflationary Hubble scale. In Ref. Q it has been shown that this 
could be achieved if the vector field is non-minimally coupled to gravity with a couplin g of the form ^RA 2 , where R 



is the Ricci scalar. This idea was employed in the vector curvaton mechanism in Ref. [lOj. Another possibility was 
to consider a varying kinetic function as discussed in Ref. Q. The first attempt to implement this idea to the vector 
curvaton mechanism is in Ref. (Tlj . 

In the above early works on the vector curvaton, attention was paid mainly on the transverse components of the 
vector field perturbations. The reason was that the vector curvaton mechanism employs the fractional perturbation 
of the density of the vector field, which is a scalar quantity, i.e. isotropic. Indeed, the mechanism does not cause any 
anisotropic expansion because, when the zero-mode of the field undergoes rapid coherent oscillations, so do the field's 
perturbations. Thus, the transverse components were deemed enough to generate ( and the longitudinal component 
was ignored. This turned out to miss a key effect introduced to the curvature perturbation from vector fields, namely 
statistical anisotropy. 

Statistical anisotropy in the curvature perturbation is a new observable and it amounts to direction dependent 
patterns in the CMB temperature perturbations as well as in the large scale structure (e.g. rows of galaxies) if it is 
intense enough. The observational limit on such a signal is rather weak as statistical anisotropy in the spectrum of £ 
is allowed up to as much as 30% [f2j . This implies that it is likely to be observed by the forthcoming observations of 
the Planck satellite mission which will detect statistical anisotropy in the spectrum if it is larger than 2% [TH, while 
there is no reason to think that this will be the ultimate limit. Another possible indication of statistical anisotropy 
in the curvature perturbation is the alignment of the quadrupole and octupole moments of the CMB (the so-called 
"Axis of Evil" [3]) and of gala xy spins [HI, [T^]. The alignment of the low CMB multipoles was shown to persist 
beyond foreground removal [l7| despite being statistically extremely unlikely with isotropic perturbations (see also 
[3). Moreover, in Ref. [l6[ it is claimed that the preferred direction of the alignment of galactic spins appears to 
coincide with the "Axis of Evil" , as one would expect if they were both due to statistical anisotropy in £. It has to 
be noted that statistical anisotropy cannot be introduced by scalar fields alone, so if its observation is confirmed, this 
would be a powerful indication that vector fields at least affect (\ 

The first study of the generation of statistical anisotropy by a vector field is in Ref. [lj|, where a vector field 
is assumed to modulate the decay of the inflaton. In a separate development [20| it was pointed out that even 
a homogeneous vector field will generate statistical anisotropy in £, because it will break the statistical isotropy 
of the scalar field perturbation generated from the vacuum fluctuation. This is further explored in Ref. [2lJ. A 
comprehensive model- independent study of statistical anisotropy can be found in Ref. [22| . The study employs the 
SN formalism [23| . extending it to include the effects of vector fields. In Ref. [22| it is shown that a vector field 
can alone generate £ through the vector curvaton mechanism only if the particle production process is approximately 
isotropic. Otherwise, the contribution of the vector curvaton field has to be subdominant and its significance amounts 
to the generation of statistical anisotropy only. As an application, the ^RA 2 vector curvaton model was revisited and 
shown to result in statistical anisotropy of order unity, which would be excessive if the vector field were to provide 
the dominant contribution to £. In such a case, the only way that excessive statistical anisotropy can be avoided 
without introducing another dominant (and statistically isotropic) source of curvature perturbations (such as scalar 
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fields) is by considering either a triad of identical vector fields each perpendicular to the other two [24| or hundreds of 
identical vector fields, randomly oriented in space. In such a case anisotropy in the spectrum or in the expansion can 
be avoided or suppressed, which means that vector fields could also act as inflatons, without generating anisotropic 
expansion. 

Usin g ve ctor fields as inflatons was first considered in Ref. [25| . More recently, vector field inflation has been studied 
in Ref. j26j with the \RA 2 model, involving hundreds of randomly oriented vector fields. For more studies of similar 
vector inflation models see Refs. [271I281. 1 In the same spirit, vector fields have been considered as candidates for dark 
energy [2^, [3(|. Furthermore, Refs. [3l| and [32| investi gate inflation or dark energy respectively, using Yang-Mills 
fields. Finally, inflation with P- forms is explored in Ref. |33j . 

Also obtained in Ref. 12211 is the contribution of vector fields in the bispectrum of the curvature perturbation (for 
the trispectrum see Ref. |34|). Using these results, in Ref. [HJ the non-Gaussianity in the spectrum is studied in 
a model- independent way, though the results are applied to the models of Refs. [10[ and [ijj. A recent study in 
Ref. [HI, extends this treatment in the case of non-Abelian vector fields, applying the results in the \RA 2 model. In 
all cases, it is shown that statistical anisotropy in the spectrum and bispectrum are correlated. This is a smoking gun 
for the contribution of a vector field to the curvature perturbation. It is likely that the Planck mission will observe 
non-zero non-Gaussianity in the CMB temperature perturbations. If this is the case then simple single-scalar-ficld 
models of inflation will be ruled out. If an angular modulation of the non-linearity parameter /nl, which quantifies 
non-Gaussianity, is also observed, then this is a strong indication of the contribution of vector fields to 

An important aspect of the mechanisms which break the conformality of vector fields and result in particle pro- 
duction of their perturbations during inflation is the appearance of instabilities. For example, in Ref. [37| the \RA 2 
model was criticised for giving rise to ghosts and other instabilities, e.g. when the modes of the perturbations exit 
the horizon, although in Ref. [22| the latter was demonstrated not to be a problem (also see Ref. [22| for a discussion 
on ghosts) . For a comprehensive analysis on ghosts and tachyon instabilities see also Ref. [38| . 

In this paper we introduce a new model of vector curvaton, which does not suffer from such instabilities (see also 
Ref. [1^]). The model consists of a massive Abelian vector field with a Maxwell type kinetic term, whose kinetic 
function is modulated during inflation similarly to the models in Ref. 0. Our model is a follow up of Ref. [llj ]. 
which did not consider statistical anisotropy. We attempt a detailed analytical investigation and confirm our findings 
numerically. As we demonstrate, our vector curvaton model gives rise to a rich phenomenology. In particular, when 
the mass of the vector field is smaller than the inflationary Hubble scale the vector curvaton can give rise to statistical 
anisotropy in both the spectrum and bispectrum, with the angular modulation dominating /nl, rendering thus the 
model, in this case, falsifiable by observations in the near future. When the mass of the vector field is larger than the 
inflationary Hubble scale though, particle production is isotropic and the vector field can alone generate the observed 
curvature perturbation. The parameter space for both regimes is shown to be substantial. 

The structure of our paper is as follows. In Sec. II we present our model and obtain the equations of motion of the 
perturbations of the vector field components in momentum space. In Sec. Ill we study particle production in order 
to obtain the desired scaling of the mass and the kinetic function of the vector field, which results in the production 
of a scale invariant superhorizon spectrum of perturbations. In Sees. IV and V we calculate in detail the power 
spectra of the perturbations in the two possible alternatives for the scaling of the kinetic function. To this end we 
obtain the evolution of the mode functions of the perturbations, through analytic approximations of high precision. 
In Sec. VI we calculate the statistical anisotropy in the spectrum and bispectrum of the produced perturbations in all 
cases. In Sec. VII we study the evolution of the zero mode of the vector boson condensate, during and after inflation. 
The evolution of the zero mode is necessary in order to determine the cosmology after inflation as well as calculate 
the curvature perturbation. In Sec. VIII we apply the vector curvaton mechanism in order to obtain constraints on 
the model parameters such as the mass of the vector field and the inflationary Hubble scale. In Sec. IX we present 
two toy-models which demonstrate how the desired scaling of the mass and the kinetic function of the vector field 
can be obtained in a natural way in the context of supergravity and superstrings. In Sec. X we summarise and 
discuss our findings. Finally, in Sec. XI we present our conclusions. Throughout the paper we use natural units such 
that c = h = fcs = 1 and Newton's gravitational constant is 8nG = rrip 2 , where nip — 2.44 x 10 18 GeV is the reduced 
Planck mass. 



1 Another way to avoid anisotropy is to consider a time-like vector field [29ll . 
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II. EQUATIONS OF MOTION 



Consider the Lagrangian density for a massive, Abelian vector field 



(1) 



d u Afj,. The above can be 



where / is the kinetic function, m is the mass and the field strength tensor is 
the Lagrangian density of a massive Abelian gauge field, in which case / is the gauge kinetic function. However, we 
need not restrict ourselves to gauge fields only. If no gauge symmetry is considered the argument in support of the 
above Maxwell type kinetic term is that it is one of the few (three) choices [38| which avoids introducing instabilities, 
such as ghosts |37ll. Also, we note here that a massive vector field which is not a gauge field is renormalizable only if 
it is Abelian [40j . 

We focus, at first, on a period of cosmic inflation, during which we assume that the contribution of the vector field 
to the energy budget of the Universe is minute and can be ignored. Thus, we take the inflationary expansion to be 
isotropic. We also assume that inflation is of quasi-de Sitter type, i.e. the Hubble parameter is H as constant. We 
will consider that, during inflation, / = fit) and m = m(t) can be functions of cosmic time t. 

Inflation is expected to homogenise the vector field. Following Ref. [ll[, we find that the temporal component of 
the homogeneous vector field has to be zero, while the spatial components satisfy the equation of motion: 



H+- \ A + —A = 



f 



f 



(2) 



where the dot denotes derivative with respect to t. From the above it is evident that the effective mass of the vector 
field is 



M = 



m 

77' 



(3) 



where we assumed that m, / > 0. 

In order to study particle production we need to perturb the vector field around the homogeneous zero mode (t) 
as: 



A^{t,x) = + SA^x) . 
Now, let us switch to momentum space by Fourier expanding the perturbations: 

d 3 k 



/d k 
( 27r )3/2 5 4«(*> fc ) ex P(* fc • x ) 



(4) 



(5) 



Then, according to Ref. j 1 11 ] the equations of motion for the spatial components of the vector field perturbations in 
momentum space are: 




5 A 1 - = 



and 



H + 



f 



2H + 2- 




m 

T 



SA n = , 



(6) 



(7) 



where a = a(t) is the scale factor of the Universe, k = \k\ and the longitudinal and transverse components are defined 
as: 



&A" = k(k : SA) and SA L ^SA-SA^ 



(8) 



To continue we need to employ the canonically normalised vector field and the physical (in contrast to comoving) 
vector field W^, whose spatial components are: 



B=^/fA and W = B/a=y r fA/a. 



(9) 



5 



Expressing Eqs. (|6|) and (JT)) in terms of the physical vector field we obtain the following equations for the spatial 
components of the physical vector field perturbations: 



and 



where 



d? + ZHd t 



1(1) _I__f_ 

2 1 f f r 



f \a 



5W 1 - = 



3H+ ( 2H + 2--1 ' 



m f) (h) + rf 



d t + - 



2\f f f 



+ 



1^2 



v 2f){ m f) { tf + ^ f 



a 



SW 11 = 0, 



/H 3 k 
——5W(t,k) exp(ifc-x). 



(10) 



(11) 



(12) 



III. PARTICLE PRODUCTION 



In this section we attempt a preliminary study of the particle production process during inflation in order to 
ascertain what kind of time dependence do f(t) and m(t) need to have in order to result in scale-invariant power 
spectra for all the components of the supcrhorizon vector field perturbations. 

To study particle production of the vector field during inflation, we first need to promote SW to a quantum operator. 
To do that we expand in creation and annihilation operators as 



/ d?k > r 

6W= J —-^3 }2[ex{k)&x(k)w x (t,k)( 



e* x (k)a{(k)w* x (t,k) ( 



ik ■ : 



(13) 



where k = k/k, with k = \k\ and A = L,R,\\ with L,R denoting the Left and Right transverse polarisations 
respectively such that 



(^W^) 2 = (5W L ) 2 + (8W n ) 
The polarisation vectors e\ can be chosen as 

ei = -L(l,i,0), e R = -L(l,-i,0), e,| =(0,0,1) 

while we have canonical quantisation with 

a x (k),a{,(k')\ = (2ir) 3 5(k-k')5 xx , . 



R\2 



(14) 



(15) 



(16) 



Since Eqs. (|10[) and (fTTj) arc linear they will be satisfied also by the corresponding mode functions w x (t,k). To 
study particle production in this theory we need to solve these equations with the appropriate initial conditions. 
After obtaining the solutions we will find the appropriate constraints on / and m, which can provide us with a scale 
invariant spectrum of superhorizon perturbations. The strategy is to begin with the transverse equation (|10p . since 
it is the simplest. Conditions obtained by requiring a flat spectrum for the transverse components of the vector field 
perturbations will hopefully simplify the longitudinal equation (fTTj) as well. 



A. The transverse components 



Let us assume that the time dependence of the kinetic function can be parametrised as 



(17) 
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where a is a real constant. We will also assume that / — > 1 at (least by) the end of inflation so that, after inflation, 
the vector field is canonically normalised. Then, according to Eq. (| 10[) the equation of motion for the transverse mode 
functions is 



wl.r + 3Hw l ,r 



1 



(a + A)(a-2)H 2 + M 2 + 



wl,r = , 



(18) 



where we used Eq. ([3]) and considered that the theory is parity invariant so that <5Wl = <$W_r, i- e - wl = wr. = Wl,r- 
In analogy to the equation of motion for the mode function of a scalar field during quasi-dc Sitter inflation, it can 
be readily deduced that a scale invariant spectrum of perturbations is attained if 



(i.e. either / oc a 2 or / oc a ) and 



-1±3 



(19) 



(20) 



where the star denotes the time when the cosmological scales exit the horizon. The latter condition simply requires 
that the physical vector field is effectively massless at that time. 2 The above condition can be obtained as follows. 
Let us define 



aM 



(21) 



If we assume that the field is effectively massless M -C k/a (i.e. r < 1) until the end of inflation we can straightfor- 
wardly calculate the power spectrum of the perturbations. First, we solve Eq. (j!8[) using the Bunch-Davis vacuum 
boundary condition 



lim wl 



^ik/aH 



2k 



(22) 



which is valid well within the horizon, where the inflationary expansion is not felt and one can use flat spacetime 
quantum field theory. The solution is then 



w L , R = SL^lJ^ e *^H) H P(k/aH) , 



where H v is the Hankel function of the first kind and v is given by 

v=-\a + l\. 

At late times (superhorizon scales) the dominant term in the above solution approaches 



lim Wl r 



ui 



-3/2 



2r(l H 



2aH 



(23) 



(24) 



(25) 



Hence, the power spectrum of the perturbations is 



Vl,r 



2^2 



lim w L r 



+ + 



4-7T 



H\ 2 ( k ^- 2 » 



[r(l - v)] 2 \2tt J \2aH 



(26) 



As evident, scale invariance is attained when v = 3/2, which, in view of Eq. (|24p. results in the values for a shown in 
Eq. (fl9|) . If a = — 1 ± 3 then the above gives 



"Pl,. 



(27) 



2 Note that this is not the same as having A M being effectively massless. In the latter case the vector field is approximately conformally 
invariant and docs not undergo particle production. However, the conformality of the massless physical vector field is broken. 
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i.e. the spectrum is scale invariant with amplitude given by the Hawking temperature for de Sitter space. 

As is shown later on, even if the r C 1 condition is violated before the end of inflation (but after the cosmological 
scales exit the horizon) , the scale dependence of the spectrum at the vicinity of the cosmological scales is not altered 
(only its amplitude is). The reason for this is that, when r> 1, the mass term in Eq. (|18p dominates the fc-dependent 
term. Consequently, the equation loses its sensitivity on scale dependence, which means that the evolution of the 
perturbations in the r> 1 regime will not affect their dependence on scale. Hence, a scale-invariant spectrum will 
remain so. We demonstrate explicitly this in Sec. IIVI Thus, Eqs. (|T9")) and (|20[) arc sufficient for the generation of a 
scale invariant spectrum of perturbations for the transverse component of our vector held. 



B. The longitudinal component 



Let us assume now that the time dependence of m can be parametrised as 

ui oc a , 



(28) 



where /3 is a real constant. Then, in view of our notation and Eq. ([lip, the equation of motion for the longitudinal 
mode function is 



"'ii 



2-a + 2^ 
1 + r 2 



Hwu 



2 - a + 2/3 
1 + r 2 



H 2 + - (1 + r 2 ) 



w\\ = . 



(29) 



Now, let us make use of Eq. (|19[) . which is necessary to obtain a scale- invariant spectrum for the transverse components. 
If we assume also that r< 1, the above equation simplifies to 



w\\+{5-a + 2/3)Hw\\ + 



-^(a-2)(2-a + 2f3)H 2 + 



Wn = . 



(30) 



Similarly to the transverse components, we can solve the above equation by using the vacuum boundary condition, 
which reads 



lim W „ = 7 ^= e lk ' aH 



2k 



(31) 



Note that, for the longitudinal component of the vector held perturbations, the vacuum condition is multiplied by the 
Lorentz boost factor 7, which takes us from the frame with k = (where there is no distinction between longitudinal 
and transverse components) to that of momentum k 7^ 0. The Lorentz boost factor is 



e _ y/(tf + M 2 _ [—£ 

M M V r 2 ' 



7= M = 



(32) 



where we considered that the effective mass of the physical vector held is given by Eq. ([3]). 

Now, assuming that r <C 1 at early times, Eq. (|30j) remains valid within the horizon, where it can be matched to 
the vacuum expression in Eq. pip with 7 ~ 1/r. By doing so we obtain 



a- 3 / 2 



"'II 



AH sin(7ri>) 



Jo 



k 



k 



(33) 



where Jo denotes Bcssel function of the first kind, and 

. 1 



v = 2 V 79 + 2 (« + 1)(2 - a + 2/3) + (2 - a + 2f3) 2 . 



(34) 



At late times (supcrhorizon scales) the dominant term in the above solution approaches 



lim w\\ 



-3/2 



7T e l f ( p +§) ( H 



r(l - 0) V H sin(7ri>) \MJ \2aH 



k 



(35) 
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where we used Eq. (|21j) . Hence, the power spectrum of the perturbations is 

167T (H\ 2 (H\ 2 ( k ^~ 2 ° 



k 3 



lim ' 



sin 2 ( 7 r£')[r(l-i>)] 2 \M J \2tt J \2aH , 
It is clear that scale invariance is attained when v — 5/2. If this is so then the above becomes 



(36) 



i.e. the spectrum is scale invariant with amplitude given by the Hawking temperature for de Sitter space, but also 
determined by the mass fraction M/H of the physical vector field. As with the case of the transverse components, 
scale invariance is retained even in the regime when r>l (see Sec. IIV|) . 

Now, using Eqs. (|TT)|) and one finds that scale invariance (i.e. z> = 5/2) requires 

/3 = -|(3±5). (38) 

We can discard one of the above values as follows. 3 According to Eqs. (TTT)) . (f2~Tj) and (j2"B")) 

r«« 1+ ^/ 2 . (39) 

This means that the value = —4 is unacceptable, because it implies that r is a decreasing function of a (for 
a = — 1 ± 3), which means that our assumption r«l cannot hold true in the subhorizon limit. Hence, we conclude 
that a scale invariant spectrum for the longitudinal component of the vector field perturbations can be attained only 
if = 1, i.e. 

m oc a . (40) 

In view of the above and Eqs. ([27]) and ([3T|) . we see that, if a = 2, then (when r < 1) we have 

V\\ = constant > Vl,r (for a = 2) , (41) 

where we considered Eq. (|20[) . On the other hand, if a = —4, then V\\ oc a~ e . Thus, even though the spectrum is 
scale invariant, its amplitude decreases in time. Scale invariance is not spoiled by the piling up of more and more 
perturbations, while they exit the horizon, because their amplitude is reduced accordingly in time, as shown in Fig. [T] 
The gradual decrease of V\\ implies that, even though {V\\/Vl.r)* ^ 1) this ratio decreases towards the end of inflation 
and may allow for the longitudinal and transverse spectra to be eventually comparable. 

IV. CASE: / oc a~ 4 & m oc a 

In this section we assume the values 

a = -4 & = 1, (42) 

i.e. / oc a -4 and M oc a 3 [cf. Eq. ((3))]. Using these we perform an analytic calculation of the power spectra for all the 
components of the supcrhorizon vector field perturbations generated by the particle production process. 

A. The transverse components 

For the choice in Eq, (|4"2l , the equation of motion for the transverse mode functions in Eq. (fT5)) becomes 

w LtR + 3Hw L ,R+(^\ (l+r 2 )w L ,R = 0. (43) 



3 The '±' sign in Eq. I|19| l is uncorrelated with the one in Eq. J38J, which means that discarding one of the values of j3 does not imply 
doing so also for one of the values of a. 



9 



■'""'\V\\ 


A \ A 


\ 


\ \ \ 

t a 







H 



a/k 



FIG. 1: Log-log plot of the superhorizon power spectrum V\\ oc a 6 (case a = —4) in terms of the physical lengthscale £ ~ a/fc 
at a given fixed time t. The spectrum is flat and it is shown by the solid horizontal lines, which depict its value at two different 
times: t\ and ti > t\ for superhorizon scales £ > H~ l . The slanted arrows show the evolution of superhorizon modes of given, 
fixed k. The figure attempts to show that, as time passes and more perturbation modes exit the horizon, their amplitude at 
horizon crossing reduces in such a way that they end up on top of the fiat spectrum at the time of exit. 



The generic solution to the above can be obtained when r 96 1. Indeed, one finds 

k 

aH 



wl,r = a 



-3/2 



w l ,r = a 3/2 



(A 

M 



C l J 3/2 [j^ ) +C2J-3/2 



for r < 1 and 
for r> 1, 



(44) 
(45) 



where Cj and Cj are integration constants. Because of Eq. (|20[) , we find that the switch-over between the two regimes 
occurs after horizon exit since 



M* k 
If < ajj 



< 1 



(46) 



and, in view of Eqs. (j39J> and (I421) . 

r oc a 4 . (47) 
Also, note that at the switch-over moment, denoted here with subscript k (because it is /c-dcpendcnt), we have 



H a k H 



«1, 



(48) 



exactly because the switch-over moment for a given mode occurs when the mode in question is already of superhorizon 
size. 

Using the above we can approximate the solutions in Eqs. (|44p and (|4"5"]) in the vicinity of the switch-over as 



-3/2 



1 



k 



3 Cl V aH 



3/2 



C2 



(aH 



3/2 



for r < 1 



and 



(49) 



10 



WL,R 



Cl 



M 
3H 



1/2 



C2 



3i/ 



1/2 



for 



> 1 



(50) 



where we used that k/a k H <C 1, M k /H <C 1 and lim J v (x) = ggrfl+^l ■ 

To approximate the complete solution we can perform a matching at the switch-over of both the mode function 
wlh and its time-derivative wl^r. In effect, this will enable us to connect the q integration constants with the c^. 
To perform the matching it is useful to consider that 



M 
3H 



oc a 



<H_ 

a 



.-1/4 



After matching we obtain 



Cl 



y/3 (2a k H 
~~T \ k 



c 2 and c 2 = — — 



2a k H 



Cl 



(51) 



(52) 



Now, to evaluate the Ci we need to match the solution in Eq. (|4"4")l with the vacuum condition in Eq. (|2"2"j) . Thus, we 
obtain 



1 / 7r in 
ci = — \ — and ci = — » / — . 

Hence, the analytic expression for the transverse mode functions can be written as 



w l ,r = a 3/2 



wl,r = 



fOT aH ~ 1 ' 
for — < 1< — 




(53) 

(54) 
(55) 



where we have used that 



(2a k HY 



16 /off 

y It 



SH 

IT 



= constant, 



(56) 



(57) 



as obtained by Eqs. ([21]) and ([5T|) . Eq. (|55|) corresponds to both the Eqs. p9|) and ([50]) after the matching, as they 
are both of the form wl,r = Ci + C2(J~ 3 , with d being constants, which are identified in the two regimes by the 
matching. Calculating the power spectrum using Eq. (|55[) . one finds the result shown in Eq. (|27[) . This demonstrates 
that the validity of this result extends beyond the r <C 1 regime of the previous section, to the region with r> 1 as 
long as M < H. 

The above analytic approximation is remarkably accurate as shown by comparison with the solution of Eq. (|43[) 
obtained numerically (see Fig. The reason is that the physical scale k/a evolves over an exponentially large range 
of values, whereas the approximation might be challenged only over a couple of orders of magnitude around the 
switch-over value k/a k = M k (i.e. r = 1). 



B. The longitudinal component 

Let us concentrate now on the longitudinal mode function. For the choice in Eq, (|42p . the equation of motion for 
the longitudinal mode function in Eq. (|2"5|) becomes 



Hw\\ 



24 
1 + r 2 



:H 2 



(1 + r 2 ) 



(en 



(58) 
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The generic solution to the above can be obtained when r 96 1 and r 76 r c (cf. Eq. (|62[) ) . Indeed, one finds 



wu = a 



wu = a 3/2 



-11/2 



C3J1 



5/2 



k 
7Jl 



5/2 



off 



^3/8 (vf^ +C ^-3/«(vf;) 



,-3/2 



M 



C3J1/2 ( 3^:] +64J-1/2 ( ^7 



3if 



for r<l, 

for l«r«r c , 

for r r c , 



where Cj, d i and q are integration constants and r c is determined by the condition 



r H 



(59) 
(60) 
(61) 

(62) 



which shows that r c 3> 1 (considering Eq. (|20[0 . The above condition corresponds to the moment when the first and 
last terms in the square brackets in Eq. (j5"8")) become comparable. 4 Note that Eq. (fBTj) looks the same as Eq. (|45l) . 
This is because Eqs. (|43|) and (|58|) become identical for r ^> r c . 

In this case we need to do two matchings of toy and w\\; at r = 1 and r = r c . Let us perform the matching at r c 
first. In the vicinity of r c Eqs. (|60|) and (|61[) can be written respectively as 



= a 



-3/2 



3 1 



3/8 



r(ll/8) \ \l 8 r 



r(5/8) 



C3 



M 
3ff 



1/2 



C4 



3H 



3/8' 



1/2' 



for r < 



for r > 



and 



(63) 
(64) 



Both the above are of the form W11 = C\ + C^a 3 , which means that they can be readily matched. Indeed, after 
matching one finds 



4 = 4 



3 /8 



7T V 3 



3/16 



r(ii/8) 



2a k H \ 



1/2 



c 4 and c'^\^- f 



3/16 



r(5/8) 



( k 



\2a k H 



1/2 



C3 



(65) 



Now we can proceed similarly to the transverse case and write the solutions in the vicinity of the r = 1 switch-over 



wu = a 



-11/2 



15 C3 \aH 



-3/2 4 /f 
7T 



C3 



M 
3F 



5/2 



1/2 



3C4 



C4 



3H 



1/2- 



for r < 1 and 



for r > 1 



(66) 
(67) 



Note that Eq. (|57| incorporates both Eqs. and (|6~4"|) . by virtue of the matching in Eq. 

As in the transverse case, to approximate the complete solution, we perform a matching at the r = 1 switch-over 
of both the mode function w\\ and its time-derivative w\\, connecting thereby the Cj integration constants with the Cj. 
After matching we obtain 



C3 



4V3 _ 4 



C3 and £4 = v3 



t 4 (T)' 



1 



rC 4 



64 



135 V 2a k H 



C3 



(68) 



4 To visualise the physical interpretation of r c consider the range [k/a,H], whose lower limit reduces in time. The effective mass of the 
physical vector field grows in time as M tx a 3 . When it becomes comparable to k/a we have r = 1 and a = a^. After this moment, M 
enters into the region in question and continues to grow approaching H. The moment where r = r c corresponds to when M crosses the 
geometric mean of the region, i.e. M = *J ' H(k/a). 
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To evaluate the Ci we need to match the solution in Eq. (1591) with the vacuum condition in Eq. (|31[) . Thus, we obtain 

1 



(':>, 



and C4 



2 K V H 2 K V H 

Hence, the analytic expression for the longitudinal mode function can be written as 



"'ii 




AY 




k 




for 




aH J 








for 


k 













> 1 



< 1 < 



~M 



(69) 

(70) 
(71) 



aH\ Z/2 /3H\ 1/2 (M 
— ImJ 1/2 [~3H 



for **<!. 
M ~ 



(72) 



In the above, Eq. (fTTj) corresponds to both Eqs. ((70)) and (f?2")) . This is because both equations, in the vicinity of the 
r = 1 switch-over, are of the form w\\ = Ca~ 3 , with C a constant (see Eqs. (|66[) and (p7| . the latter also incorporating 
Eq. (|60p in the limit r > 1). This corresponds to the growing mode of Eq. (|70|) and the decaying mode of Eq. I|72[). 
which is dominant near a ^ ak- It can also be shown that this mode remains dominant until M ~ H, i.e. beyond the 
r ~ r c switch-over. Calculating the power spectrum using Eq. (jTTj) , one finds the result shown in Eq. (j3"7| . As with 
the transverse case, this demonstrates that the validity of this result extends beyond the rCl regime of the previous 
section, to the region with r ^> r c as long as M <gi H . As was the case of the transverse modes, the above analytic 
approximation is remarkably accurate which is demonstrated by comparison with the solution of Eq. (|58p obtained 
numerically (see Fig. [5]). 



C. The power spectra 



In the previous section we found that the transverse and the longitudinal mode functions evolve in a similar manner. 
Indeed, the functions cease oscillating after horizon exit and until M ~ H. During this period they satisfy an equation 
of the form (cf. Eqs. ([50) and ([67])) 



= «-3/»4/f 
7T 



CA 



M 
3H 



1/2 



CB 
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1/2- 



for M, -<H, 
a 



(73) 



where the c, are appropriate constants. The above expression is of the form w\ = C \ + C^a -3 , with C,; constants. 

However, there is a crucial difference between the transverse and longitudinal modes. In the case of Wl r, ca arL d eg 
in Eq. (|73[) arc identified with c\ and ii respectively, which are shown in Eq. (|52p . Using these values it can be readily 
confirmed that the dominant term in Eq. (|73[) in the regime M,k/a <C H is the "growing" mode: wl,r = constant. 
In contrast, in the case of w», ca and cb in Eq. ()73|) are identified with £3 and £4 respectively, which are shown in 
Eq. (|68|) . Using these values it can be readily confirmed that the dominant term in Eq. (|73|) in the regime M, k/a <C 
is the decaying mode: twii cx a -3 (by virtue of the condition M/H <C 1 <?C aH/k). 

This difference in scaling is reflected, of course, on the power spectra, as already shown in Eqs. (j27|) and (|37|) . which, 
as argued above, are valid in the regime where M, k/a H . Thus, the typical value for the vector field perturbations 
in the regime in question scales as 



SWl,r 
where SW\ = \5W\\ and 



H 

2^ 



constant 



and 



SWu 



3H H 
~M2^ 



cx a 



(74) 



-(Vl+Tr) =V l ,r , 



(75) 



because Vl = Vr as the theory is parity invariant. 
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FIG. 2: Log-log plot of the evolution of the real part of the longitudinal (upper) and transverse (lower) mode functions wa 
with A = ||, L, R of the vector field perturbations at a given, fixed k, in terms of the physical momentum scale k/aH (weighted 
by the Hubble scale H) in the case when a = —4 (/ oc a -4 ). The plots have been normalised with respect to the amplitude of 
the transverse mode functions at horizon crossing (k/aH = 1). The solid lines depict the numerical solutions of the equations 
of motion in Eqs. (|43p and (|58|l . while the dashed lines depict that corresponding analytic approximations in Eqs. (|54[) - (|56p 
and ([70} - (|72|) respectively. The precision of the approximation is remarkable as the difference can hardly be seen. Notice 
that there are clearly three regimes for the evolution of the mode functions: First, when subhorizon, they undergo oscillations 
until horizon crossing (k/aH ~ 1), when they enter a power-law regime, which eventually is terminated by another phase of 
oscillations, when the vector field becomes heavy (M > 3H). The longitudinal mode function, when superhorizon, scales as 
W|l oc a~ 3 before oscillations, in agreement with Eq. (|71(l . In contrast, the transverse mode function remains constant with 
respect to a, as suggested by Eq. (|55|l . Both mode functions oscillate with amplitude ||wa|| oc a~ 3 when the vector field becomes 
heavy, cf. Eq. ((82)) . 
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Now, let us find out how the perturbations behave when M ~ H. Eqs. (|56|) and (|72j) can be written as 

~R ( fc \ 3/2 i , , , 1 



where we defined 



2 \i—\jj) w L , R = -^=J 1/2 {z) + -x z ^J„ l/2 {z) and (76) 

/ k\ 3/2 1 1 

2 \J-[h) w \\ = ^'M/a^-^J-^W, (77) 

z = — oc a J and x = — < 1 . (78) 
6H an 

When z > 1, for the amplitudes of the oscillating Bessel functions we have ||Jj,(2:)|| « ||J_ v (z)||. Using this and also 
that i < 1 for superhorizon perturbations, we obtain 



|ii>l,r|| ~ ||w|| || for z > 1 . (79) 

Indeed, when z> 1 the mode functions approach 

i { H\ i/2 smz 1 fH\ 3/2 cosz 

WL ' R= vm\k) — and w " = ~vifUJ ~ forz>>1 ' (80) 

where we used that lim J^^) = y^j cos[z — ^-(^ — |)]. Since z>l, the frequency of oscillations is very large 

compared to the expansion rate H. This means that it makes sense to use the average values of the power spectra 
over many oscillations. In view of Eq. (|80[) . we find 

1_ /£ x 



^+ = ^ll=^2(^) ■ (81) 



Thus, the typical value for the vector field perturbations in this regime is 

where we used 5W\ = s/V\, where A = +, || and we have SWl = SWr = 8W + . We have found that, for M > H, the 
mode functions of the transverse and longitudinal components oscillate rapidly with the same amplitude, so that the 
components of the typical perturbations of the vector field are equal. 



D. Deviations from scale invariance 



So far we have discussed the case of exact scale invariance for both the transverse and longitudinal spectra. A small 
deviation from scale invariance is favoured at the moment by the observations. Such a deviation can be achieved by 
perturbing the scalings of / and m. Indeed, introducing the perturbations 

a = -4(l + e/) and /3 = 1 + e m , (83) 

one can repeat the calculations to find the scale invariance of the spectra. Parameterising this scale invariance in the 
usual manner, i.e. 

V x oc fc"*-\ (84) 
with A =||, +, wc have found the following expressions for the spectral indexes: 

n\\ -l = -(e/ + 2e m ) and n+ - 1 = --(18e/ + e m ) . (85) 

We see that, for positive e^'s, both the spectra are red as required by the observations. The spectral indexes are 
identified if e rn sa 2ef. In this case n\ ~ 0.96 is obtained with Cf ~ 8 x 10 -3 . Deviation from scale invariance can 
be also produced by considering that inflation is not exactly de Sitter, i.e. e = —H/H 2 ^ 0. According to the usual 
curvaton result, this would introduce a red tilt in the spectral indexes of magnitude 8n\ = —2e. Thus, e w 0.02 may 
be enough to satisfy the observations, for (a, 0) = (—4, 1). 
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V. CASE: / oc a 2 & m oc a 

Here we study the other possibility which allows scale invariant spectra for all the components of the perturbations 
of our vector field. Hence, in this section we assume the values 

a = 2 & /3 = 1, (86) 

thus making / oc a 2 and M = constant. The condition that the physical vector field is effectively massless at the time 
of horizon exit in Eq. (f20| suggests that M/H <C 1 at all times when the scaling above holds. Using this condition 
and scaling we can calculate the power spectra for all the components of the superhorizon vector field perturbations 
generated by the particle production process. 

A. The transverse components 

For the choice of scaling above, the equation of motion for the transverse mode functions in Eq. (fTS")) becomes 

fc^ 2 



w l ,r + 3Hwl,r+ I -J {l + r^WLM^O, (87) 

which has the same form as the equation of motion for the a = —4 case. However, because M is constant we may 
obtain an exact solution to the equation of motion above 



-3/2 ' n 

1 w 



where 



1 / fM\ 2 3 



PS 2V 9 " 4 U "2- (89) 



As before, we have used the Bunch-Davis vacuum boundary condition in Eq. (|22|) to obtain the constants of integration. 
Taking the limit of this solution at late times (superhorizon scales) we can calculate the power spectrum 



H ,2 



n=(^], 00) 



v ~ | — | (M/H) , then we get a slight deviation from scale invariance V+ oc k n + 1 where 



2 [M\ 

n+-l = - -EF = 2r lA- (91) 



which is in agreement with Eq. (|27p and we considered v ~ | . However if we consider the next level of perturbation 

3 _ 1 

2 3 

2 

3 \H { 

As evident, this deviation from scale invariance gives rise to a blue tilt on the spectrum. Alternatively we may perturb 
the scaling slightly as 

a = 2(l + e/ ) and f3 = 1 + e m , (92) 
in analogy with Eq. (|83j) . If e/, m > M/H then the scale dependence becomes 



n+ -l = --e/, (93) 



which allows a red spectrum if e/ < 0. 
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FIG. 3: Log-log plot of the evolution of the real part of the longitudinal (upper) and transverse (lower) mode functions w\ 
with A = ||, L, R of the vector field perturbations at a given, fixed k, in terms of the physical momentum scale k/aH (weighted 
by the Hubble scale H) in the case when a = 2 (/ oc a 2 ). The plots have been normalised with respect to the amplitude of the 
transverse mode function at horizon crossing (k/aH = 1). The solid lines depict the numerical solutions of the equations of 
motion in Eqs. (|87|) and Q94p . while the dashed lines depict the corresponding analytic solution in Eq. (|88p and the approximation 
in Eqs. (J9TJ) and (|98[) respectively. The precision of the approximation is remarkable as the difference can hardly be seen. Notice 
that there are clearly two regimes for the evolution of the mode functions: First, when subhorizon, they undergo oscillations 
until horizon crossing (k/aH ~ 1), when they enter a power-law regime. Both mode functions, when superhorizon, remain 
constant with respect to a in agreement, e.g. with Eq. (|98|) for wu. Because M <C H, we have w\\ S> wl,r- 
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B. The longitudinal component 



Let us concentrate now on the longitudinal mode function. For the choice of scaling in Eq. (|86|) . the equation of 
motion for the longitudinal mode function in Eq. (|29l) becomes 



3 + 



Hwu 



(1+r 2 ) 



wu = . 



1 + r 2 

Using the same technique as in the previous case we find exact solutions in two regimes 



(94) 



-5/2 



"'H 



for r <C 1 
c 5 + CqO,~ 3 for r 1 



and 



(95) 
(96) 



where a, £i are integration constants. Again we have used M/H = constant < 1. As in the previous case, using the 
vacuum boundary condition in Eq. pip , we can find the values of 05^- Then, matching the two solutions and their 
first derivatives at k/aH = M/H (i.e. r = 1) we can find 65^ in the regime r> 1. Thus, the solutions become 



H 2 



/2k Mk 



2 (M 
31+ T 5 [h 




aH 



fOT aH ~ * : 



(97) 
(98) 



As in the previous section, the above analytic approximation is remarkably accurate, as demonstrated in Fig. |3l 
Using Eq. (|98p , we can calculate the power spectrum at arbitrary late times 



(99) 



which is in agreement with Eq. (|3"?| . Since in this case, we have M/H = constant <C 1, we see that the longitudinal 
power spectrum is constant, in contrast to the a = —4 case. Also, we see that V\\ ^> V+. 

Once again if we introduce perturbations to the scaling as in Eq. (|92p we can find a slight deviation in the scale 
dependence 



'II 



1 = 2( £/ - e rn ) . 



(100) 



VI. STATISTICAL ANISOTROPY AND NON-GAUSSIANITY 

The theory studied in this paper has two clear advantages. First we can obtain a completely isotropic perturbation 
spectrum for the vector field, which has previously never been achieved. As we discuss below, this means that we 
may consider vector fields as dominating the total energy density of the Universe when the curvature perturbation 
is formed. The second advantage is that we can also account for a small amount of statistical anisotropy in the 
curvature perturbation spectrum depending on when inflation ends, again by considering the vector field alone. We 
also demonstrate this in what follows. Finally, statistical anisotropy can also be present in a correlated manner in the 
bispectrum as well, which characterises the non-Gaussian features of the CMB temperature perturbations. In view 
of the forthcoming observations of the recently launched Planck satellite mission this is a particularly promising and 
timely result. 

Let us begin with the generic treatment before focusing on cases of particular interest. As in Ref. [ia |. keeping only 
the leading term in the anisotropy, the power spectrum may be parameterised in the following fashion 



P c (k) =Pf 



(k) 



1 + 9 



(d-ie) 



(101) 



18 



where P 1 ^ is the isotropic part of the power spectrum, k = k/k, and d is a unit vector in some chosen direction. In 
the above the parameter g quantifies the statistical anisotropy in the spectrum. In Ref. [l2j . it has been shown that, 
currently, the data sugg est that g < 0.3 or so. 

As shown in Ref. |22J , when the anisotropy in the spectrum is due to the contribution of a vector field, then at tree 
level (first order) we have 

d = N A , (102) 
where Na = Na/N a , Na = \Na\ and the components of Na arc defined as 

dN 

with N being the number of remaining e-folds of inflation. The vector Na quantifies (to first order) the contribution 
to the curvature perturbation £ from the vector field only; as defined in the SN formalism (23[. In this formalism the 
magnitude of anisotropy is characterised by |22j 

J^i^GT ( 104 ) 



■JJ1SO 



and the isotropic part of the spectrum is 



Vf° = N 2 ^ + N A V+ = Np><, ( 1 + f ^ ) 



where we have defined 



and 



(106) 



dN 

N^-. (107) 

In the above <f> is a scalar field which also contributes to the curvature perturbation. Its contribution is necessary 
if the anisotropy due to the vector field is excessive. In this case the vector field can only generate a subdominant 
contribution to Q, with the dominant contribution due to some other source, e.g. the scalar field </>. If, however, 
the anisotropy due to the vector field is within the observational bounds then one can dispense with the scalar field. 
Thus, if this is the case, we can assume that the scalar field modulation of N is negligible so that — > 0. Therefore, 
dispensing with the scalar field contribution is equivalent to £ — > oo. 

Now, from Eqs. (| 105[) and (|106[) we can write the anisotropy parameter in Eq. (|101[) as 



9 = t: 



£P 4 



(108) 



Hence, we can calculate the level of statistical anisotropy in the power spectrum. 

How do the above translate in the case of the vector curvaton mechanism? As calculated in Ref. [22], the tree- level 
contribution to the curvature perturbation spectrum from the vector curvaton field is 



4Q\ 



V + + (V\\-V + )(W-kf 



where W = W /W and W = \W\. In the above &a is defined as 



n A 



PA 

p 



(109) 



(110) 



where pa is the density of the vector field. As in the case of the scalar curvaton paradigm, Eq. (|109[) should be 
evaluated at the time of decay of the curvaton field. 
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As we have shown, in the case when a = —4 and if M > H by the end of inflation, V+ ~ V\\ and statistical 
anisotropy in the spectrum can be very small; within the observational bounds. This means that the vector field can 
alone give rise to the observed curvature perturbation. Thus, if this is the case, we can dispense with the scalar field 
<f) and we can write P c (fe) = V CA {k). Then, from Eqs. ([TUT]) . (TTU5|) and (fTP9"|) we find 



7^° = N%V+ with N A = ~± 



and also g = ( r P\\/'P+)_— 1 which agrees with Eq. (|108p in the limit £ — > oo. 

As shown in Ref. [22], apart from statistical anisotropy in the spectrum of the curvature perturbations a vector 
field can give rise to statistical anisotropy in the bispectrum, which characterises the non-Gaussian features of the 
cosmological perturbations. Non-Gaussianity is an observable of great importance as it can be a powerful discriminator 
between mechanisms for the generation of the curvature perturbation. It was shown recently in Ref. (35| that the 
statistical anisotropy in the spectrum and bispectrum are correlated if they are generated by a vector field. Such 
correlation can be a signature prediction of the vector curvaton scenario. If it is observed it would be direct evidence 
for a vector field contribution to the curvature perturbation. 

Non-Gaussianity expresses a non-vanishing 3-point correlator, or equivalently a non-zero bispectrum. The bispec- 
trum is a function of three fei,2,3 vectors which may be chosen arbitrarily. However we will focus on the two standard 
configurations; the equilateral and squeezed (or local) configurations where the magnitudes satisfy k\ = lt2 — k% and 
k\ ~ fca 3> &3 respectively. It was shown in Ref. [35j that the expressions for /nl, the non-linearity parameter charac- 
terising non-Gaussianity, for a scale invariant power spectrum in the equilateral and local configurations respectively 
are given by 



6 f c quil ,2^2 3 (l + |g 2 ) + [p+|(p 2 -2g 2 )]iyl 

5 /nl ~ 5 l + 2h A 



and 



6 .local _ C 2p2 3 1+pWj 



where W± = \W±\, with W± being the projection of the unit vector W to the plane defined by the three fci,2,3 
vectors. 

We have also defined 

V\\ - V+ V- 
V = — 1 and q=^~, (114) 

with 

P- = ^(V l -Vr) = 0, (115) 

i.e. g = since our theory is parity invariant. 

In the following, to calculate /nl, we assume that V<p is due to a light scalar field, in which case 

"♦'(I) 2 - <116 > 

The curvature perturbation in Eq. (|10ip and the vector field contribution to it in Eq. (|109[) should be calculated 
at the time of the decay of the vector field, as in the scalar curvaton scenario. In contrast, the spectra of the vector 
field perturbations Vx are calculated either at the onset of the vector field oscillations or at the end of the scaling of 
/(a) and m(a); whichever occurs latest, regardless of whether this moment is during or after inflation. Therefore, at 
that time, / = 1 and M = m, where rh = constant is the final (vacuum) value of the vector field mass. 

As we have shown, the evolution of the components of the vector field perturbations goes through the following 
stages. When subhorizon, the perturbations are oscillating (starting of as quantum fluctuations). After horizon 
crossing, they follow a power law evolution of the form 5W\ = C\ + C2a~ 3 , where Ci are constants. In the a = —4 
case the effective mass M is increasing with time, which may allow the the vector field and its perturbations to begin 
oscillations (when M ~ H) before the end of the scaling behaviour of / and m, i.e. before the end of inflation. In 
contrast, in the a = 2 case, since M = constant <C H , the perturbations remain in the power-law regime until the end 
of scaling. This suggests that, in this case, rh < H . Thus, when a = 2, the oscillations commence only after inflation, 
when H(t) decreases enough such that rh ~ H(t). In view of the above we can now calculate the statistical anisotropy 
in the spectrum and bispectrum of our model. 
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A. The power-law regime 

The curvature perturbation is formed in this period of evolution if rh ^ -ff*, where denotes the Hubble scale 



during inflation. Then, from Eqs. and (|37J) we have V\\ = 

our model gives 



V+, which means that V\\ > V+. Therefore, 



3H t 



and 



- 1 



(117) 



where £ is to be calculated at the decay of the vector field, which occurs after the onset of the vector field oscillations. 
The latter, in this case, has to occur after the end of inflation. 

In this regime we also observe that V+ = V^. Substituting Eq. (|117[) into Eqs. (|112p and (|113[) (with M — > rh) we 
can now get expressions for the non-linearity parameter /nl in terms of the anisotropy parameter g. We find 



G 



/•cquil 

/nl 



3 








H 




2Q A 







(M^) 2 -l 



20a 



/3if, 



V rh 




3iJ* 



2 V H^l 



2 V WX 



(118) 



and 



u rlocal 

5 /nl 




In analogy with Eq. (|101[) we can write 



/nl = (l + G ■ Wl) 



(119) 



(120) 



where Q quantifies the statistical anisotropy in the bispectrum. It is evident that the direction of the statistical 
anisotropy in the bispectrum is correlated with the one in the spectrum (cf. Eq. (|109p ). since they are both determined 
by the direction of the unit vector W exactly as in Ref. (35|. At the moment, observations do not provide any 
information about the value of Q . 



Highly anisotropic and subdominant limit 



Suppose that, at the end of scaling, V\\ 
p = (3H*/rh) 2 > 1. Thus, Eq. fTUS]) gives 



3> V+. This is equivalent to (H^/rh) 2 3> 1. In this limit we find 



9 = 



3iL 



(121) 



Were we to dispense with the scalar field, i.e. if the vector field contribution to C were dominant, then £ 3> 1 and the 
above would give g = (3H*/rh) 2 3> 1, which clearly violates the observational constraints. Therefore, the contribution 
of the vector field perturbations to the curvature perturbation must be subdominant to that of the scalar field. Thus, 
we must have Na <C N^, which means ^ < 1. Similarly, the vector field contribution to the total energy density of 



the Universe has to be subdominant as well fl A "C 1. 
when (H*/rh) 2 > 1 wc find 



Hence, Eq. (|110[) suggests CIa — > j^a- Using this, for the case 



G 



•cquil 
NL 



2g 2 / rh 
04 \3H~* 



1 / 3ff» 
8 V rh 



wl 



(122) 
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and 

'3£L 



6 f local „ V 

5 /NL VI a V3iJ* 



wl 



(123) 



From the above we see that, in this case, /nl oc g, i.e. statistical anisotropy in the spectrum of the curvature 
perturbation would intensify the non-Gaussianity as is the case in Ref. [35|. Furthermore, by comparison with 
Eq. (|120p . we see that, for both equilateral and local configurations, for the isotropic part we have 

i^-^im) = i? (124) 

which can, in principle, be substantial even if ^ <§C 1 because VIa can be very small. Notice however, that, in both 
configurations, Q 3> 1, by virtue of the condition rh < H*. This means that, in this case, the statistical anisotropy in 
the bispcctrum is dominant and would be readily detected if there is a confirmed detection of a non-zero /nl- Turning 
this around, if non-Gaussianity is not observed to be predominantly anisotropic this would rule out the rh < regime 
of this model. 

From Eqs. (|122[) and (|123|) we readily obtain that the amplitudes of the modulated /nl in both configurations are 



i\\f^\\=2^-(^L\ and -||/^ uil H = -— (125) 



Hence, in general 1 1 1 1 < II/nl"'!!- This means that the observed upper bound /nl < O(10 2 ) [4l[ should be applied 
to the equilateral amplitude. Hence we obtain that g 2 /VLa < 10 2 , which means that /^£ < 0(1), since rh < H* in 
this regime. 

Before moving on, it is important to stress that the above results are valid for a = — 1 ± 3, i.e. for both / oc a~ 4 
and / cx a 2 , as the condition rh 2 <C H 2 is possible in both cases. In fact, this is the only possibility for the case when 
a = 2. 



2. Almost isotropic and dominant limit 

In the case when a — —4, M is growing with time and the bound in Eq. (|20| can be satisfied even with rh not 
smaller than . Therefore, in this case we can investigate the possibility that rh ~ , which corresponds to the 
edge of the power-law regime (which is inaccessible in the a = 2 case) . 

Considering this regime we find V\\ ~ "P+ = V$. Thus, in this case we find an almost isotropic curvature perturba- 
tion, which allows us to dispense with the scalar field and assume that the curvature perturbation is dominated by 
the vector field contribution. Hence, we can take Na 3> N^, i.e. £ 3> 1. Then, from Eqs. (|108|) . (|114[) and (|117[) we 
find 

3HA 2 ST , , 

m J V+ 

where 5V = V\\ — V+. Thus, if the fractional difference of the spectra is not excessive, the vector field might generate 
statistical anisotropy in the CMB within the observational bounds. Using the above, it is straightforward to show 
that 

| /NL = _l_( 1+9 ^), (127) 

in both the local and equilateral configurations. By comparison with Eq. f|l 20[) we see that, in this case, Q ~ g. Hence, 
statistical anisotropy is of the same magnitude in both the spectrum and bispectrum of the curvature perturbations. 
Therefore, observational constraints on g in Ref. (l2T ] suggest that /nl may feature an angular modulation at a level 
as large as 30% or so. 

In the isotropic and dominant limit when m — > 3H* we find that g = and 

/nl = -4- . (128) 

As expected, this is equivalent to the scalar curvaton scenario @- Thus, in this case, substantial non-Gaussianity can 
be generated if the vector field decays before it dominates the Universe, with PIa -C 1. 
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B. The late-time oscillations regime 



As we have shown in Sec. IIV CI in the case when a = —4, it is possible that M grows much larger than H before 
the end of scaling. This is indeed so if rh ^> . Then, after M ~ H, the vector field perturbations undergo rapid 
oscillations, during which, V\\ = V+ as shown in Eq. (|8Tj) . This means that g = and p = (cf. Eqs. (|108|) and (|114|) 
respectively), which results in an isotropic power spectrum. In this regime, from Eqs. (|81[) and (I116p . we also find 
that 



1 / 3H N 

^»<+=2 M) V *- 



Therefore Eqs. {TT2]) and ([TT5|| reduce to 

6 



NL 



2fl A 



( rh 



\3H* 



(129) 



(130) 



which is valid in both the local and equilateral configurations. 

Since the vector field perturbations are isotropic we have no need of the scalar field contribution. As discussed 
previously, we can dispense with the scalar field in the limit £ — > oo. In this limit Eq. (|130|) reduces to the scalar 
curvaton expression in Eq. (|128p . as expected. However, if there is indeed a contribution to the curvature perturbation 
from a scalar field, this may affect the value of /nl even if £ 3> 1. The reason is easily understood from Eq. (|129p . 
which suggests that, at the end of scaling V\\.+ <C V^, since M-»m> iJ*. 5 Thus, if (m/3i?*) 2 > £/2 3> 1, we find 
from Eqs. (fTT2"j) and (|TTg|) 



, ,2 5 /P+ 
/nl = 4 77— ^~ 



40, 



£ ^3iJ + 
2 



< 



40a 



(131) 



which is smaller than the scalar curvaton result. 



C. One-loop corrections 

In the above we considered only the tree level contribution to the spectrum and bispectrum of the vector field. In 
principle, one-loop corrections may also contribute significantly. Their contribution in the case of vector fields has 
been studied in Ref. [42j , where it is shown that the one- loop corrections to spectrum and bispectrum dominate only 
if 

N 2 AA V A > Nl (132) 

where V A = (2V+ + V\\) and N AA = \N AA \, with 

d 2 N 
dW l dW, 

As shown in Ref. [22|, for the vector curvaton we have 



A A = aTTTTTFiF ■ (133) 



— and N A a = 2 — ^ 

3 W AA W 

Using the above we can recast Eq. (|132[) as 



N A=ol£ and N aa = 2^. (134) 



6W~ifp^>\w, (135) 

which clearly violates our perturbative approach. Hence, we conclude that in the vector curvaton case the one-loop 
correction has to be subdominant. 



5 This reflects the fact that the amplitude of the oscillating vector field perturbations becomes exponentially suppressed during inflation 
because it decreases as a~ i before the end of scaling, cf. Eq. I I82II . 
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VII. EVOLUTION OF THE ZERO MODE 



As is evident from Eq. (|109[) , in order to calculate the curvature perturbation associated with the vector field one 
needs to study also the evolution of the homogeneous zero mode W. Combining Eqs. © and © and using Eq. (|17l) . 
we obtain 



W + 3HW 

where we also considered Eq. 



1 - -a) H--(a + 4)(q - 2)H 2 + M 2 
2 / 4 



w = 0. 



(136) 



A. During inflation 

As we have shown, to obtain a scale invariant spectrum for the transverse components of the vector field pertur- 
bations we require /(a) to scale according to Eq. (|19[) . i.e. a = — 1 ± 3. Using this and considering de Sitter inflation 
(with H « 0) the above becomes 6 



W + 3HW + M 2 W = 0. 



(137) 



We show below that, when M«iJ (true at early times when a = —4; always true when a = 2), the solution of the 
above is well approximated by 



W ~ & + C 2 a 



(138) 



where Ci are constants. The dominant term to the solution of Eq. (|138|) is determined by the initial conditions. A 
natural choice of initial conditions for the vector field zero-mode can be based on energy equipartition grounds. As is 
demonstrated in what follows, if energy equipartition is assumed at the onset of inflation, the dominant term turns 
out to be the decaying mode, W oc a~ 3 when a = —4 and the "growing" mode W = constant, when a = 2. 

To apply energy equipartition in the initial conditions we need to consider the energy-momentum tensor for this 
theory, which, from Eq. |T]), is given by [ll| 



( \g^F pa FP° - F^fA + m 2 (a„A v - \g^A p A" 



If we assume that the homogenised vector field lies along the £ 3 -direction, we can write the above as pd| 

T£ = diag(p A , -p_L, -p±, +p ± ) , 

where 

PA = Pkin + Va , PJ_= Pkin - Va , 

with 



Pkir 



V A 



--fF,uF^ = -a~ 2 fA 2 = - 
4 p 2 2 



W 



1 ~ \ a 1 



m 2 A,,A^ = -a- 2 m 2 A 2 



1 



M 2 W 2 , 



where A = \A\, we used Eqs. © and p7|) . and we assumed a negative signature for the metric. 
Energy equipartition, therefore, corresponds to 

(pidn)o - (Va)o , 

where the subscript '0' indicates the values at some initial time, e.g. near the onset of inflation. 



(139) 

(140) 
(141) 

(142) 
(143) 

(144) 



The equation of motion for the zero mode of the canonically normalised field is: B + HB + (-2H 2 + M 2 )B = 0, which agrees with the 
findings in Ref. 
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1. Case: f oc a 



In this case M oc a 3 and the solution to Eq. (|137[) is 



W 



l M 



I M 
C 2 eos( — 



(145) 



When M > H the above shows that the amplitude of the oscillating zero mode is decreasing as ||W|| oc a 3 . In the 
opposite regime, when M <C H the solution above is well approximated by Eq. (|1 38[) with C\ = C^a^ 3 Mo/SH , where 
we considered that a~ 3 M = Mo = constant. Using this, the constants C2 and C3 in Eq. f|145[) can be expressed 
in terms of initial values of the field amplitude Wq and its velocity Wq (in field space): 

Wo , (w Q +3HW ] 



C 2 = 



3H 



and C3 



Mo 



i . 



(146) 



Assuming initial equipartition of energy we can relate Wq with Wo- From Eqs. (|142[) and f|143[) . setting a = —4, 
we readily obtain 



Then, using Eq. ()144[) . we get 



Pkin 



^(W + 3HW) 2 and V A = ^M 2 W 2 . 



(147) 



Wo^ W (-3H±M 



(148) 



Substituting this relation into Eq. (| 146)) we find that the evolution of the vector field W in Eq. f|145[) takes the simple 
form: 



W = W< 



a 



I M 7T 

V2 cos ± - 

' 3H 4 



(149) 



Note that this equation is valid for any value of M. However, we can see that when M <gi H the zero mode of the 
vector field is decreasing as W oc a -3 , but when M 3> H it oscillates rapidly with a decreasing amplitude proportional 
to a~ 3 . On this basis we can assume that the typical value of the zero mode during inflation always scales as 



Woe a- 3 . 



(150) 



With the assumption of initial equipartition of energy for the vector field at the onset of inflation we can calculate 
the kinetic and potential energy densities. 7 Inserting Eq. (|149[) and its derivative into Eqs. (|142[) and (|143|) we find 



Pku 



W M sin ( — ± - 
1 3H 4 



and Va 



W M cos I — ± - 
1 3H 4 



(151) 



Hence, the total energy density is constant 



p A = M 2 W^. 



(152) 



Because this relation is independent of the vector field mass M it is valid in both regimes: when M <§; H and W 
follows a power law evolution, and when M 3> H and W oscillates. This is valid as long as /(a) and m(a) are varying 
with time. 



7 By "potential" we refer to the energy density stored in the mass-term Va = — ^m 2 A^A 11 . 
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2. Case: f oc a 2 

In this case, M = constant, which means that the solution of Eq. (|137[) is 



W = a- 3 / 2 C 1 aVi Z (W + c 2 a-yn z (¥y 



(153) 



Since in this case M <C H , the above solution is always well approximated by Eq. (|138|) and there is no oscillating 
regime. 

Now, Eqs. (fT42|) and (fT43j) take the form 



Pkin = \w 2 and V A = \m 2 W 2 . 



-H 2 C 2 a- 6 - 



Combining Eqs. f|138[) and (|154p . we find 

9 

Thus, assuming energy equipartition at the onset of inflation (cf. Eq, (|144[1 ) gives 

2 



1 + 
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— » 1 => d c ±—0,0*0 



M 



where we used that Mq = M <C H . Inserting the above into Eq. (|138[) we find 



W = On 3 C 2 



constant ~ Wo , 



(154) 



(155) 



(156) 



(157) 



because, after the onset of inflation, (ao/a) 3 -Cl^ 3H/Mo- 

Therefore, we have found that W remains constant. Since M ~ constant, this means that Va also remains constant. 
On the other hand, Eq. (|155l) suggests that pkin oc a~ 6 . Thus, since we assumed energy equipartition at the onset of 
inflation, we find that, during inflation, pkin <C Va- Hence, 



Pa ~ Va — M 2 W \ 



(158) 



where M = constant = Mo- This result is the same as in the case / oc a~ 4 in Eq. (|152[) . 

We should stress here that, according to the above, for a = — 1 ± 3 the typical value of the zero mode scales as 
W oc a*-"/ 2 ^ 1 , which means that the zero mode of the comoving vector field A = aW/\f] remains constant. 



B. After the end of scaling 



The simplest choice would be to assume that the scaling of /(a) and m(a) during inflation is terminated at the end 
of inflation. This would imply that / and m are modulated by a degree of freedom which varies during inflation, e.g. 
the inflaton field. We explore this possibility in Sec. IIXI Here we briefly comment on the possibility of allowing the 
scaling to end before the end of inflation. 

At the end of scaling / — > 1 and m — > rh = constant. Then, the zero-mode equation of motion can be obtained by 
setting a = in Eq. (|136p . During de Sitter inflation this gives 



W + 3HW + {2H 2 + rh 2 )W = . 



(159) 



which is solved by 



W = a-*' 2 Idia^^ + C^a-Vi^ 



(160) 



The end of scaling occurs before the onset of the oscillations if rh <C H . 
-l 



In this case the 



where C- arc constants. 

above suggests that W = C' x ar l + C' 2 a~ 2 , with C[ constants. This is not of the same form with Eqs. (|138[) or (|153[) . 
Therefore the transition at the end of scaling is not smooth but it is indeed felt by the zero mode if this occurs before 
the onset of the oscillations. 
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The end of scaling occurs after the onset of the oscillations if m > H (which is possible only in the case with a = —4). 
In this case Eq. fl!60|) suggests that the amplitude of the oscillations decreases as ||W|| cx a~ 3 / 2 , which agrees with the 
findings of Ref. This is not in agreement with Eq. (|145p . which suggests that, during the oscillations, ||W|| oc a~ 3 . 
Thus, even if the end of scaling takes place after the onset of the oscillations, the transition modifies the evolution of 
the zero-mode. 

The vector field perturbations mimic the behaviour of the zero mode if the end of scaling occurs in the r» 1 
regime, i.e. when rh 3> k/a. For such scales, after the end of scaling, the equations of motion for the transverse and 
longitudinal components, Eqs. (|10|) and (fTTj) . become identical and of the same form as Eq. (|159|) . Thus, the mode 
functions satisfy the equation 

w A + 3Hw x + (2H 2 + m 2 )w x = , (161) 
which means that the typical value of the vector field perturbations is 



8W\ = \[V\ ocw\ = a 3/2 ciaViMw) 2 + c 2 a W^ST 2 



(162) 

where Ci are constants. Therefore, similarly to the above, the solution is not similar to Eqs. (f5"5j) and (fTTj) for the 
transverse and longitudinal components respectively if rh <C H. This is also true if rh > H, i.e. when the end of 
scaling occurs after the onset of the oscillations, for which ||dWx|| oc a -3 / 2 ; in contrast to ||<5W,\|| cx a -3 before the 
end of scaling cf. Eq. (|82|) . Note here that the spectrum for the scales that exit the horizon after the end of scaling 
is not scale invariant. 

From the above we see that, if the end of the scaling occurs before the end of inflation the evolution of the zero- 
mode and the vector field perturbations is non-trivially affected, which significantly complicates the treatment. Since 
attributing the scaling of f(a) and m(a) to some degree of freedom which varies during inflation but not afterward 
(e.g. the inflaton field, see Sec. IIX[) is much more physically motivated than the alternative considered above, in the 
following, we assume that this is indeed the case and the end of scaling occurs at the end of inflation. 



C. After inflation 



At the end of inflation we assume that the scaling of / and m has ended and we have 

/ = 1 and m = rh . (163) 

Hence, Eqs. (|152[) and (|158|) no longer apply. The evolution of pa is determined as follows. 

As mentioned already, after the end of scaling, a = and M = rh. Then, Eqs. (|142[) and (| 143[) become 

P^ a = \{W + HWf and V A = \rh 2 W 2 . (164) 

The behaviour of pkin arid Va depends on whether the vector field is light or not. To sec this let us calculate the 
evolution of the field after inflation. With the conditions in Eq. (1163j) the physical vector field of Eq. (j9)) is W—A/a, 
while Eq. (|136|) becomes 



W + 3HW +IH + 2H 



where the Hubble parameter after inflation decreases as H(t) 
of the Universe. Solving Eq. (|165j) we find 



m 

2 



W = 0, (165) 
3(i-t-w)t ' w ith w = p/p being the barotropic parameter 



W 

W + HW 



CiJdirht) + C 2 J-d{rht) 



and 



CiJd-i {hit) - C 2 J\-d {rht) 



(166) 
(167) 



where d = g^ji 3 ^ ■ One can easily see that the vector field behaves differently if it is light, rht <C 1, or heavy, rht ^ 1. 
Let us first see what happens if the vector field is light. Then, Eqs. ()166[) and (|167[) can be approximated as 
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Although the solution has one decaying and one growing mode, it might happen that the decaying mode stays larger 
than the growing mode. To check this we calculate constants C\ and C2 by matching the above equations to the 
values Wend and W cn d at the end of inflation (denoted by 'end'). Thus, we find that 



W 



a 

3W + 1 V Gend 
-2 



i(3™-l) 



W + HW = H* 



^end 



end 
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and 



(170) 
(171) 



where if* is the inflationary Hubble scale. Plugging these solutions into Eq. (|164p (and using that a 3< - 1+w > cx t 2 ) we 
obtain 



Va 

Pkin 



(3w + l) 2 



t 



_] (_) ~(m^<l: 



which implies that the total energy density of the light vector field is 



Pa — Pki] 



end 



Q-cnd 



pA oc a 



(172) 



(173) 



Therefore, we see that the energy density of the light vector field scales as that of rclativistic particles. This is in 
striking difference to the scalar field case, in which when the field is light its density remains constant even after 
inflation. 

On the other hand, if the vector field is heavy, rht ^> 1, the Bessel functions in Eqs. (|166[) and (|167p are oscillating 
and the latter can be recast as 



W + HW 
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7T 



Ci cos ( rht 



l + 2d 



7r I + Co cos f rht — 
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7T 



C\ sin rht — 



l + 2d 



tt) + C2 sin rht 



1 - 2d 



1 - 2d 



and 



(174) 
(175) 



As can be seen above, when the vector field is heavy, it oscillates with a very high frequency and with amplitude 
decreasing as t _1 /( 1 +"') cx a~ 3 / 2 . This was already demonstrated in Ref. [l[. When we calculate the energy density 
of the oscillating vector field from Eqs. (|141[) and (|164[) we find 



PA 



-rht !+« 



C\ +Cl + 2(71(72 cos (dir) 



Pa oc a 



(176) 



Thus, we see that the energy density of the heavy vector field scales as that of non-relativistic matter. Furthermore, 
if calculating the pressure from Eq. (|14ip . we find that it oscillates with a high frequency: 



1 A 2_ 

= —mt 1+" 

7T 



C n sin (2mt - dn) + (7f sin (2mt + dir) + 2C X C 2 sin (2mt) 



(177) 



Therefore, we have found that, on average, the oscillating vector field behaves as pressureless isotropic matter (see 
Eq. (|140p ) . which is in agreement with Ref. [l|. Hence, the massive vector field can dominate the Universe without 
generating excessive large scale anisotropy. This is crucial for the vector curvaton mechanism because, to produce 
the curvature perturbation, the field must dominate (or nearly dominate) the Universe without inducing anisotropic 
expansion. 



VIII. CURVATON PHYSICS 



In this section we calculate constraints for our vector curvaton model assuming that the scaling behaviour of /(a) 
and m{a) ends when inflation is terminated. This implies that the scaling is controlled by some degree of freedom 
which varies during inflation, e.g. the inflaton field. A specific example of this kind is discussed in Sec. IIXI 
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A. Basics 



According to the curvaton scenario [2( the total curvature perturbation can be calculated as the sum of individual 
curvature perturbations from the constituent components of the universe multiplied by the appropriate weighting 
factor. In our scenario this is written as follows 



C = (1 - A )Crad + ti A ( A , (178) 

where &a is defined in Eq. (|110j) . As in the scalar curvaton paradigm, the above is to be evaluated at the time of 
decay of the curvaton field. 

As was discussed in Sec. E2 if m ^> at the end of inflation, then the vector field perturbation spectrum is 
isotropic and may generate the total curvature perturbation in the Universe without violating observational bounds 
on the statistical anisotropy of the curvature perturbations. If this is the case, we can assume that £ rac i = 0. On 
the other hand, when to <C H*, the amplitude of the spectrum of the longitudinal component of the vector field 
perturbations is substantially larger than the one of the transverse perturbations. Hence, the curvature perturbation 
due to the vector field is excessively anisotropic. To avoid conflict with observational bounds (see Ref. [12]), the 
contribution of the vector field to the curvature perturbation has to remain subdominant. Therefore, for this scenario, 
we have to consider £ rac i ^ and the curvature perturbation already present in the radiation dominated universe must 
dominate the one produced by the vector curvaton field. 

For definiteness let us assume that inflation is driven by some inflaton field, which after inflation oscillates around 
its VEV until reheating, when it decays into relativistic particles. The vector field must be subdominant during this 
time. But after reheating, the Universe is radiation dominated with the energy density decreasing as /9 ra d oc a~ 4 . If 
the vector field at this epoch is heavy and therefore undergoes rapid oscillations, its relative energy density increases, 
PA/pr&d oc a, as can be seen from Eq. (|176[) . When the field becomes dominant (or nearly dominant) it can imprint 
its perturbation spectrum onto the Universe. 

The contribution to the curvature perturbation by the vector field is calculated as follows. On the spatially flat 
slicing of spacetime we can write for each component of the content of the Universe 

C n = -H^-, (179) 

Pn 

where n represents different components of the cosmic fluid. Using the continuity equation p n — — 3H (pn ~f~.Pn), one 
can recast the above equation for the vector field as 



, Sp A 

Qa = - — 

OPA 



„ 2 \\SW\\ 
dec ^3 \\W\\ 



dec 



2 SW 

3 W 



(180) 

end 



where we considered that the decay of the vector field (labelled by 'dec') occurs after inflation and after the onset of 
its oscillations so that it is pressureless, as shown in Eq. (|177p . In the last relation we took into account that, after 
inflation, the equations of motion are the same for the zero mode and for the superhorizon perturbations of the field. 

In Sec. ITVCl it was shown that the typical value of the field perturbation is 8W ~ (3H*/M)(H t ,/2n). If M < H* 
this is because the longitudinal component is dominant over the transverse ones (see Eq. ([74)) '). If AI ^> H*, then the 
transverse and longitudinal components are oscillating with the same amplitudes (see Eq. (|82p). 

For this reason, at the end of inflation, we can write 

6W end ~^^d$, (181) 
to Zir m 

where we have taken M = rh and / = 1 at the end of inflation. IU on d can be found from Eq. (|152[) by using 
(p A )cnd ^ WqMq ~ W cnd m (see Eqs. (fl52j) and ([158]) ). Thus, 



Wcnd „ _ (182) 

TO 



Hence, from Eq. (|180[) we calculate the curvature perturbation of the vector field 



Ca-^- 1 / 2 — , (183) 
mp 
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where f2 e nd = (pa/ p) end is the density parameter of the vector field at the end of inflation, p on d is the total energy 
density dominated by the inflaton field, and we have used the Friedmann equation: 'SmpH^ = p C nd- Since the vector 
field must be subdominant during inflation we have Q cn d "C 1. Since we also require (a ~ (SW/W) cn d < 1 for our 
perturbative approach to be valid, we obtain the following range for O cn( j: 

(—) <O end <l. (184) 
\nip J 

Eqs. (|183|) and (|184[) are valid in both a = — 1 ± 3 cases. The only difference is that, in the / oc a 2 case, statistically 
isotropic curvature perturbations cannot be generated. Hence, only the considerations for statistically anisotropic 
perturbations in Sec. lVIIIB2l arc relevant. 8 



B. The parameter space 

Here we calculate the parameter space for the model following the method of Ref. [ll[. First, we note that at the end 
of inflation the inflaton field starts oscillating and id/-1. Therefore the Hubble parameter decreases as H(t) ~ 
In general, the inflaton potential is approximately quadratic around its VEV. Thus, the coherently oscillating inflaton 
field corresponds to a collection of massive particles (inflatons) whose energy density decreases as a -3 . When the 
Hubble parameter falls bellow the inflaton decay rate T, the inflaton particles decay into much lighter relativistic 
particles reheating the Universe. After reheating, the Universe becomes radiation dominated with energy density 
scaling as p ra d oc a~ 4 . 

On the other hand, the evolution of the energy density of the vector field, depends on its mass m. As discussed in 
Sec. IVII1 if rh <C the energy density scales as pa oc a -4 until the vector field becomes heavy and starts oscillating. 
If m 3> i?*, however, the vector field has already started oscillating during inflation and pa oc a~ 3 . 

To avoid causing an excessive anisotropic expansion period the vector field must be oscillating before it dominates 
the Universe and decays. This requirement implies that 

r, m > T A , H dom , (185) 

where Ta is the decay rate of the vector field and ifdom is the value of the Hubble parameter when the vector field 
dominates the Universe if it has not decayed already. Working as in Ref. [ll| , we can estimate -ffdom as 

/ r art i ^ r 7ti *i ^ / ^ 

tfdo m ~^ri/Vin 1;- minjl; — j. . (186) 
Similarly, if the vector field decays before it dominates, the density parameter just before the decay is given by 

1/2 r - 2/3 r . ^ -1/6 

m \ . m 1 



fi dec ~ Q end I — 1 minjl;— j minjl;— | . (187) 

where fidec = (^A)dcc- Combining the last two equations and using Eq. (|183p we can express the inflationary Hubble 
scale as 

rj ( - i -1/3 ( t n 1/12 / it 1 \ I/ 4 

H * X /2 /- . „;„ J 1. m I ™;„Ji. m l /^maxjrASiidomP 



-j-^-Cx^ll;-! m „|l;-J ^ ' j| j ■ (18S) 

The bound on the inflationary scale can be obtained by considering that the decay rate of the vector field is 
T^i ~ h 2 m, where h is the coupling to the decay products. Then we can write max {r^; i/dom} ^ h 2 rh. Furthermore, 
we must consider the possibility of thermal evaporation of the vector field condensate during the radiation dominated 
phase. If this were to occur, all the memory of the superhorizon perturbation spectrum would be erased. The 
thcrmalisation rate is determined by the scattering rate of the massive boson field with the thermal bath and is given 



The lower bound in Eq. 118411 guarantees that 8W/W <g 1 throughout inflation. The reason is the following. As discussed before 

Eq. CEEQ, during inflation 5W/W ~ H^/MW. Now, for / cc a 1±3 and m oc a we have M = m/y/J ex a 2 ' ' . Also, from Eqs. (j 150 ^ 

and ||157|I we obtain W oc a5' _1±1 '. Thus, we see that, in all cases considered MW = constant, which means that SW/W = constant 
during inflation. Therefore, (SW/W) cn( ± <C 1 is sufficient to guarantee the validity of our perturbative approach throughout inflation. 
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by r sc ~ h 4 T. Requiring that the condensate does not thermalise before it decays, r sc < Ta, in Ref. [lj| it was shown 
that h must satisfy 

/ - \ 1/4 

£ h < . (189) 

rap \mp J 

The lower bound in the above is due to decay through gravitational interactions, while the upper bound becomes 
irrelevant if the vector field dominates the Universe before it decays (then h < 1 is sufficient). This is because, in the 
latter case, the energy density of the thermal bath is exponentially smaller than pa and the vector field condensate 
does not evaporate. 

From Eq. (|188p one can see that the parameter space is maximised if the Universe undergoes prompt reheating 
after inflation, i.e. if T — > H*. To find the parameter space we investigate two separate cases: when fa 3> H„ and 
when fa <^ . 



1. Statistically isotropic perturbations 

This possibility can be realised only in the case when a ~ —4. As mentioned before, if the mass of the vector field 
at the end of inflation is larger than the Hubble parameter, fa > H*, then the field has started oscillating already 
during inflation. In this case the amplitudes of the longitudinal and transverse perturbations are equal and therefore 
the curvature perturbations induced by the vector field are statistically isotropic. We can assume, in this case, that 
the vector field alone is responsible for the total curvature perturbation in the Universe without the need to invoke 
additional perturbations from other fields. Thus, we can set £ ra d = in Eq. (JT75J) and write 

C ~ ttdccCt • (190) 
Using this and the lower bound on h we find from Eq. p89[) the lower bound for the inflationary Hubble parameter 



Tj / /• \ 4/5 / 5 \ 3/5 

H* ( (, \ ( ra x 



mp ~ V vHdec / \mp 



(191) 



where we have taken into account that the parameter space is maximised when the universe undergoes prompt 
reheating, i.e. r — >• H*. From this expression it is clear that the lowest bound is attained when the vector field 
dominates the Universe before its decay, f2dcc — ^ 1, and when the oscillations of the vector field commence at the very 
end of inflation, i.e. fa — >• H*. With these values we find the bounds 

H, > 10° GcV K 1/4 > 10 14 GeV , (192) 

where U* 1 ^ 4 denotes the inflationary energy scale and we used that £ ps 5 x 10~ 5 from the observations of the Cosmic 
Background Explorer. 

In view of the above, we can obtain a lower bound for the decay rate of the vector field. Indeed, using Eqs. (|189[) 
and (ggg) we find 

771 ^ TJ 

T A > —5- > — f > 10~ 9 GcV . (193) 
ra p rap 

From the above we find that the temperature of the Universe after the decay of the vector field is 
Idee ~ V m P^A > 10 4 GeV, which is comfortably higher than the temperature at Big Bang Nucleosynthesis (BBN) 
Tbbn ~ 1 MeV (i.e. the decay occurs much earlier than BBN), and also higher than the electroweak phase transition, 
i.e. the decay precedes possible electroweak baryogencsis processes. 

Since fa > H*, Eq. (|192p corresponds to a lower bound on fa. An upper bound on fa can be obtained as follows. 
Because, fa > > T, Eq. (|187p becomes 

r~r~ 

^doc ~O endA /— . (194) 

V 1 A 

From Eq. (| 189[) we have Ta ^ fa 3 /mp. Combining this with the above we obtain 

fa* < ( ^) Tm 2 P . (195) 



dec 
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Now, when a — —4 we have Moco 3 during inflation. Since the end of scaling occurs when inflation is terminated, 
for a < a„ nr \ we can write 



M ~ e 6 "°°*H* , (196) 



where we considered that the field begins oscillating when M ~ -ff* and A^ osc is the number of remaining e-folds of 
inflation when the oscillations begin. Inserting the above into Eq. (|195|) we find 



Jv osc ^ JV ogc _ 



9 V QdecH, 

where, in the last inequality, we used that fi en d < 1 and r < H*. Now, considering that rh > if*, Eq. (|19ip gives 

> c?. (198) 

rap 

Hence, combining Eqs. (|197|) and (|198|) we obtain 

fl£r<~lnC = 4.4. ( 199 ) 

Thus, in view of Eq. (j!96|) . we obtain the bound m < e 3Ar °™ H*. which results in the following parameter space for 
rh: 

1 < rh/H., < 10 6 , (200) 

where wc used Eq. (|199[) . The above range is reduced if the decay of the curvaton occurs more efficiently than through 
gravitational couplings, i.e. if h > rh/mp. Nevertheless, we see that the parameter space in which the vector field 
undergoes isotropic particle production and can alone account for the curvature perturbation, is not small but may 
well be exponentially large. Indeed, repeating the above calculation with V a ~ m (i.e. h ~ 1) it is easy to find that 



Hence, using that ri C nd < 1 and T < we obtain 



Odec / V H 



(201) 



AC* = -|lnO dec < 3.1 => 1 < m/H, < 10 4 , (202) 

where we used fidec ^ 10 -2 . This is because, in the case considered, /nl is given by Eq. (|128|) . so a smaller rL ec 
would violate the current observational bounds on the non-Gaussianity in the CMB temperature perturbations [41| . 

Still, it seems that, to obtain an exponentially large parameter space for rh, we need pA not to be too much 
smaller than V* during inflation and also inflationary reheating to be efficient. In the case of gravitational decay 
(r^4 ~ rh 3 /mp) Eq. (|197j) has a weak dependence on both il cn d and T: moc (Q^^T) 1 ^ 3 , which means that the 
allowed range of values for rh remains large even when f2 on d and T are substantially reduced. This is not necessarily 
so when Ta ~ h 2 rh, with h 3> rh/mp. Indeed, in this case it can be easily shown that rh oc h~ 2 n^ nd T. Therefore, if T 
is very small it may eliminate the available range for rh. Fortunately, the decay coupling h can counteract this effect 
without being too small. 



2. Statistically anisotropic perturbations 

If the vector field is not responsible for the total curvature perturbation in the Universe, the parameter space is more 
relaxed. In this case, the vector field may start oscillating after inflation and hence its mass is rh -C However, this 
means that the curvature perturbation due to the vector field is strongly statistically anisotropic. For this reason we 
can no longer set C ra d to zero because the curvature perturbation present in the radiation dominated Universe must 
be dominant. In other words, the parameter £ defined in Eq. (|106p needs to be very small, £ <C 1. 
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From the Eq. (|74[) it is clear that if inflation ends when M = rh <C i?* , the longitudinal power spectrum is much 
greater than the transverse ones, V\\ 3> V+. Thus, from Eq. (| 104[) we find 



-fsr-. ( 203 ) 

where we also considered Eq. f| 105[) with (<1. Since the anisotropic contribution to the spectrum is subdominant 
[l2T ] , to first order we can write 

— -- — ^ ~ . (204) 

Hence, 

iViW| ~ K 2 - (205) 
Now, from the SN formalism at tree level we have [22| 9 

C = N,p5(j) + N A SW. (206) 
Comparing the above with Eq. (|178[) we can equate the contributions to £ from the vector field. Thus, we find 

2V A <W = Cl A U ■ (207) 

Using that 



SW 2 Y,x 8W% 2V+ + V\\ 



(208) 



<5W| <W| 
and combining Eqs. (|205[) and (|207l) we obtain 

C~.g- 1/2 n A a. (209) 

Inserting this expression into Eq. ()188[) and considering again that the lowest decay rate of the vector field is through 
gravitational decay, max {r^; -ffdom} > rh 3 /mp we find 

tt / z-2 \ 3 / 4 / - \ 5/8 / r x -3/8 r . % 1/8 

ti* / .9 C \ ( m \ ' ( I \ . | . m ' 



> TT- min 1;^ . (210) 

m P V"dec/ \ m P/ \ m -P/ I r J 

The above suggests that the lower bound on is minimised for prompt reheating with T — > H«. Also, from obser- 
vations we know that the statistically anisotropic contribution to the curvature perturbation must be subdominant. 
Thus, the vector field should not dominate the Universe before its decay. Hence, using r — > H * and Qa < 1 we obtain 



^>VHcJ±. (211) 
nip V m.p 

From this expression it is clear that the parameter space for is maximised for the lowest mass value. The 
minimum mass of the vector field can be estimated from the requirement that the field decays before BBN. Because 
the lowest decay rate is the gravitational decay, this condition reads rh? /mp > Tg BN /mp, with Tbbn ~ 1 MeV, which 
corresponds to rh > 10 4 GeV. Using this, we find that the parameter space for the vector curvaton model with the 
statistically anisotropic curvature perturbations is 

H*> ^ 10 7 GeV & U* 1/4 > g^HO 13 GeV , (212) 



Technically, here we should have NaSW —> Na-SW but the dot-product is dominated by one term, that of the longitudinal component, 
since 5W ~ SWu as shown in Eq. 12081 . 
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i.e. it is somewhat relaxed compared to the statistically isotropic case (cf. Eq. ([192])) depending on the magnitude 
of the statistical anisotropy in the spectrum, for which g < 0.3 [HJ. This result is valid for both a = — 1 ± 3 cases. 
From the above it is evident that there is ample parameter space for the mass of the vector field 

10 TcV < rh < H* . (213) 

We can readily use the above to briefly discuss an indicative example. Suppose that rh ~ 10 -2 ii*. Then, from 
Eq. (|211j) we obtain H* > 10~ 2 g( 2 mp. Assume now that statistical anisotropy in the power spectrum is observed 
with g ~ 0.01. Thus, we obtain > 10~ 13 mp ~ 100 TeV. Let us assume also that non-Gaussianity is observed with 
predominant angular modulation that peaks in the equilateral configuration with amplitude II/niI"'!! ~ 10- Then, 
according to Eq. (|125[) . we have g 2 /fldec ~ 10, i.e. ild oc ~ 10~ 5 . Also, it is easy to check that ft cn d ~ fldec^/F A /in 
comfortably satisfies the bound in Eq. (|184p if Ta > T BBN /mp. Finally, from the above we can estimate the amplitude 
of the non-Gaussianity in the local configuration ||/^ al || ~ 10~ 3 [cf. Eq. ([TZ5]) ] and also /jj s £ - 10~ 7 [cf. Eq. ipg])]. 

C. Constraints from isocurvature perturbations 

As discussed, when rh <C H* , particle production is strongly anisotropic and the curvaton has to be subdominant 
when it decays. The fact that the curvaton decays while subdominant allows the possibility to generate a sizable 
isocurvature perturbation, which needs to comply with observational constraints. Of course, if the curvaton decays 
early enough into relativistic particles which join the preexisting thermal bath then there is a possibility that no 
isocurvature perturbation is generated. For this to be so we require all the components of the late time Universe to 
be relativistic and in thermal equilibrium so that they can be produced by the decay products of the vector curvaton 
field. Hence, dark matter needs to be thermal. For WIMP dark matter we require curvaton decay to occur before 
the breaking of electroweak unification (i.e. the temperature at curvaton decay should be T dec > 1 TeV) so that the 
decay products can produce WIMPs. In this case baryons are also generated by the curvaton decay products and so 
are neutrinos. However, if the vector curvaton decay occurs later or if the dark matter is not thermal (e.g. axions) 
then an isocurvature perturbation can be generated. An estimate of the maximum isocurvature perturbation is then 
given by 

2 2 

~ p[(Cbe£ + ^dccCO — Cbcf] = rfidecCi > (214) 
before ° b 

where ildec <C 1, (Sp/p) before ~ §Cbef corresponds to the perturbation which preexists curvaton decay, while 
(5p/p)after ~ f (Cbcf + ^ccCa) [cf Eq. (|178|) ] corresponds to the perturbation after curvaton decay which includes 
the contribution to the perturbation due to the curvaton. The above assumes that one of the constituents of the Uni- 
verse content (e.g. dark matter) does not receive any contribution from the curvaton decay products while another 
such constituent receives 100% contribution, hence this is the maximum isocurvature perturbation (Realistically, S is 
model dependent and it is determined by the branching ratio of the vector curvaton decay to the various constituents 
of the Universe content). 

Observational constraints require that an uncorrelated isocurvature perturbation cannot exceed 11% of the adiabatic 
mode [43| . This means that 

fidecCl < 0.1C => 9 < 0.01 , (215) 

where we considered also Eq. (|209[) . Thus, even the maximum possible isocurvature perturbation produced by our 
model allows the generation of statistical anisotropy at the level of a few percent. 

It is important to note here that the isocurvature perturbation is uncorrelated to the adiabatic one in contrast to 
the scalar curvaton model. This is because, in our model, when rh -C -H*, C * s n °t generated by the vector curvaton 
field but it preexists curvaton decay. The vector curvaton contribution to C is negligible (i.e. C ~ Cbef) and its effect 
amounts to generating statistical anisotropy only. In the scalar curvaton model, however, it is the curvaton field that 
generates C (i-e. Cbef is negligible), which means that, when curvaton decay occurs before domination, the isocurvature 
perturbation is fully (anti)corrclated with the one of the curvaton. 

IX. SCALAR FIELDS AS MODULATORS 

Throughout this paper we have taken the modulation of the kinetic function / and the mass of the vector field for 
granted and assumed that it is due to some degree of freedom which varies during inflation. The most natural choice 



S - S ± 



Sp 
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for such a degree of freedom is, of course, the inflaton field itself but other choices are also possible. In this section 
we briefly explore a couple of such possibilities, inspired by beyond the standard model theories such as supergravity 
or superstrings. 



A. The inflaton as modulator 



In string theory the modulation of parameters such as masses or kinetic functions is due to so-called moduli fields. 
The moduli arc scalar fields which parametrise the size and shape of the extra dimensions. In that sense they are not 
fundamental scalar fields, but appear so from the view-point of the 4-dimensional observer. Typically, the dependence 
of masses and couplings on canonically normalised (i.e. with canonical kinetic terms) moduli fields is exponential. 

Consider the following kinetic function 

/(0)cce-^ /mp , (216) 
where fx > is a real constant. Comparing this with Eq. (fl~T|) we see that 

m«o«e-> p , (217) 
where we considered also Eq. (|40|) . From the above we readily obtain 

4> = -pHm P , (218) 

which suggests that 

(^) kin = \j> 2 = i M 2 p (219) 

and we have used the flat Friedmann equation p = 3(Hmp) 2 . 

Now, suppose that <fi is also driving quasi-de Sitter inflation, with ff ~ constant. This requires 
p = ~ V{4>) > {p<p)kin, which demands 



/' 



< V6 . (220) 



To drive quasi-de Sitter inflation a scalar field needs to follow a slow-roll attractor solution, in which the acceleration 
term <p in its Klein-Gordon equation of motion is negligible. Thus, the field equation for the inflaton is the well known 
slow-roll equation 

3H(j>~-V' ((/>), (221) 

where the prime denotes derivative with respect to <fi. 10 Combining the above with Eq. (|218p and taking 
V((j>) ~ 3(Hm P ) 2 we obtain 

V{4>) ~ Voe^ /mp oc a"^ 2 , (222) 

where Vq is a density scale and we used Eq. (|217[) . Thus, we find that the inflaton is characterised by an exponential 
potential, which is reasonable for a modulus field. The above also suggests that, for quasi-de Sitter inflation, when 
V{4>) ~ p0 = p ~ constant, we need to have 

fi < 1 , (223) 



10 Because pa <SC V* we can ignore the backreaction of the vector field to the scalar field dynamics. Indeed, if the modulation of / and m 
are due to a single scalar field <j> then the Klein-Gordon equation of motion of the latter obtains a source term of the form 

fir i j-f ff 

— = --f'F^F^+mm'A^ = T ( Q p kin - 2V A ) ~ — ~rPA , 

dtp A: (j) (f> 

where C is given in Eq. JT} and wc considered / <x m a and m oc a. If p\ is small enough then the above can be negligible compared 
with the —V' source term in the equation of motion. 
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which is somewhat stronger than the bound in Eq. (|220[) . Combining Eqs. (|221[) and (|222l) we get 

A0 = -m P In ^1 - ^/j, 2 AN\ , (224) 

where AiV = HAt is the number of the elapsing e-folds during which the inflaton varies by A</>. Notice that the above 
is consistent with Eq. (|218|) but it is valid only when /i 2 <C 2/AiV. 11 The violation of this condition is equivalent 
with the violation of the "robustness" assumption for H (i.c H ~ constant). That is, it corresponds to a significant 
variation in H, i.e. —AH ~ H 2 At. Technically, this can occur even under slow-roll conditions if At is large enough. 
Since a "robust" H is necessary for the generation of a scale invariant spectrum of perturbations, it should last at 
least as much as it takes for the entire range of the observable cosmological scales to exit the horizon. This range 
corresponds to about AN ~ 9 e-folds, which induces the bound 

V2 

C<y, (225) 

which is tighter than the bounds in Eqs. (|220[) and (|223[) . 

Let us attempt to obtain an estimate of \i. The slow roll parameters for the model in Eq. (1222[) arc 

e = - H /H 2 ^\ mp (pj 2 = \f, (226) 

W = m l\ = M 2 - (227) 

If the inflaton is the source of the dominant contribution to the curvature perturbation of the Universe then, the 
spectral index of is 

n s -l = 2? ?0 - 6e = -[? = -2e . (228) 

If, on the other hand, the dominant contribution to the curvature perturbation of the Universe is due to the vector 
curvaton (only possible in the a = —4 case with rh > H t ) then we have 

n s - 1 = 2r) A - 2e w -2e , (229) 

where we considered that t\a is negligible. Indeed, in almost all cases studied in the previous sections we had r\A = 
when e/, e m =0, the reason being that the power spectra for the vector field perturbations are exactly scale invariant, 
when inflation is considered to be do Sitter (i.c. H = constant), which means that any deviations from scale invariance 
would be due to f ^ 1 only. The only exception was the transverse component in the case when a = 2, where the 
T)a parameter (defined in Eq. (|91[0 may not be negligible. However, in the following, we assume that M <C H is 
sufficiently small to ignore t\a even in this case. 

From Eqs. (|228[) and (|229p we see that, regardless of whether the dominant contribution to the curvature pertur- 
bation is from the inflaton or from the vector curvaton field, the spectral index is 

n s = 1 - 2e . (230) 

Comparing the above with the observed value [4l| n s — 0.960 ± 0.014 (for negligible tensors) we find 

l-i = 0.200 i 0.036 , (231) 

which satisfies the bound in Eq. (|225[) . In view of Eq. (|224[) , the above value guarantees that H remains "robust" for 
about A7V ~ 50, which comfortably encompasses the cosmological scales. 

Note that, for the model in Eq. (|222l) . the slow-roll parameters in Eqs. (|226[) and (|227[) remain constant and smaller 
than unity. Hence, inflation never ends. One can remedy this by either invoking some kind of hybrid mechanism to 
remove the inflaton from the slow-roll trajectory and send it to the true vacuum or by modifying the model such that 
slow- roll naturally breaks down eventually. 12 



It turns out that Eq. I|218|l remains always valid even if Eq. {224) is not. 

For example, in the toy-model with inflaton scalar potential V(4>) = Vo[cosh(ii<f>/mp) — 1] the slow-roll parameters are 
r\ = i/i 2 +e and e= i/i 2 coth 2 (/i0/2mp). The model approximates Eq. (12221) when iji>mp/(i but allows inflation to end at 

<t>cnd — —rap In ( ) ~ \/2mp, where £(</> on d) = 1 and we considered fi 1 in the last equality. 
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B. Higgsed vector curvaton 



Another possibility is that the vector field is Higgsed, which means that m oc ip, where ip is a Higgs field. Hence, 
from Eq. (|40l) wc require ip oc a during inflation. Suppose that ip is rolling down a hilltop potential of the form 



V(<p) = Hop - \m 2 y + ■ ■ ■ , (232) 

where the ellipsis corresponds to terms of higher order which stabilise the potential but are negligible during inflation 
and Vtop is a constant density scale. The equation of motion of tp is 13 

(p + 3H(p - m%<p ~ , (233) 

whose solution has the following growing mode 



^a-JI^V'+iW^, (234) 
Therefore, the requirement that p oc a is achieved if the effective mass of the Higgs field during inflation is 

m v = 2H . (235) 

It turns out that such a value for the effective mass of the scalar field during inflation is quite reasonable in the 
context of supergravity theories. Indeed, in Ref. J3] it has been demonstrated that Kahler corrections to the scalar 
potential arc expected to give a contribution of order H to the masses of scalar fields. Hence, these fields would be 
fast-rolling during inflation down the potential slopes. 

Moreover, in this case, / is the gauge kinetic function, which, in supergravity theories, is a holomorphic function 
of the scalar fields of the theory. Hence, it is natural to expect that / = f(<p) and the rolling tp would modulate the 
kinetic function as well as the mass. Indeed, if m v — 2H, then to satisfy Eq. (|17|) one simply needs f(<p) oc p> a with 
a = -1±3. 

Of course, / may also depend on other fast- rolling scalar fields / = /(</>i, 02, • ■ • 4> n )- If we assume that / cx J\™ 
with on = 0(1) then 

Since Kahler corrections result in mf ~ H 2 , we find that /// ~ H as required. We would still need to tune on and 
m 2 such that /// = (—1 ± 3)iJ, but it is evident that the required values can be naturally attained in the context of 
supergravity, as also discussed in Ref. (llj . 

Having said that, since / is the gauge kinetic function in this case, we have / ~ 1/e 2 , where e oc a~ a l 2 is the gauge 
coupling. This means that, because / — > 1 at the end of inflation, the vector field becomes strongly coupled during 
inflation if a > 0. Thus, the physically motivated case is only a = —4, where e oc a 2 and the vector field is always 
weakly coupled. Fortunately, it is this case which has the richest phenomenology, as we have shown. 



X. SUMMARY OF RESULTS 



The model studied in this paper, albeit simple, has been shown to have a rich phenomenology, without suffering 
from any instabilities. By studying particle production in this model we have found that a scale invariant spectrum of 
vector perturbations for the transverse components can be obtained if a = — 1 ± 3 and the field is light at horizon exit. 
We also found that, to get a scale invariant spectrum of vector field perturbations for the longitudinal component we 
additionally require that m oc a. We have assumed that this scaling continues throughout inflation as would be the 
case if both / and m are modulated by some degree of freedom which varies during inflation, e.g. the inflaton field. 
We also assumed that inflation is quasi-de Sitter. 

In all cases we have solved the equations of motion for the mode functions of the components of the vector field 
perturbations both numerically and through analytic approximations, normalising them appropriately at the vacuum. 



As in the previous subsection we ignore the backreaction from the vector field. 
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A comparison between the two has shown that our analytic expressions approximate the numerical solutions with 
very high precision. We have found that, after horizon exit, the mode functions w\ of the components of the vector 
field perturbations cease oscillating and follow a power-law evolution of the form 

w x =C A + C B a- 3 , (237) 

where Ca,Cr are constants determined by the vacuum initial conditions. Depending on the case considered, the 
dominant term in the above is either the constant or the decaying mode. 
The effective mass of the vector field during inflation is [cf. Eq. ([3J] 

M=-^oca 1 ~( Q / 2) . (238) 

As mentioned, scale invariance for the transverse spectra requires a = — 1 ± 3 plus the vector field needs to be light 
when cosmological scales exit the horizon M* <C -ff* , where -ff* ~ constant is the inflationary Hubble scale. In the 
case when a = 2 (i.e. / cx a 2 ) we see that M = constant so that rh = M* <C ff*, i.e the vector field remains always 
light until the end of inflation. In this case, the power-law regime for the evolution of the mode functions of the vector 
field perturbations is valid until inflation terminates, while the dominant term in Eq. (|237l) is the constant one for all 
wx (see Fig. [3]). 

In contrast, when a = —4 (i.e. / cx a -4 ) we have M cx a 3 . This allows the possibility to have M* <C ff* fh, i.e. 
while the field is light when the cosmological scales exit the horizon it can become heavy before the end of inflation. 
If this is the case then the power-law regime is terminated before the end of inflation when the vector field becomes 
heavy and begins coherent oscillations. Hence, when a = —4 we could have either rh < -ff* or rh > ff*. Studying the 
power-law regime we have shown that the dominant term in Eq. (|237p for the mode functions is the constant term 
for the transverse components but not for the longitudinal one, for which the dominant term is the decaying mode 
w\\ cx a -3 (see Fig. [2]). If m < then inflation ends before the power- law regime is concluded. In this case, the 
results are the same as in the a = 2 case. However, if rh > ff* the power-law regime is terminated before the end of 
inflation. 

Consider first that rh < -ff*, which is possible for a = — 1± 3. As mentioned, in this case the power-law regime 
continues until the end of inflation. If this is so we have shown that the power spectra for the transverse and 
longitudinal components of the vector field superhorizon perturbations are given by 

ff*\ 2 , „ /3ff*\Vff* N 



v+= \^) and ^< = UM te) ■ (239) 

where V\\ is the spectrum of the longitudinal component, V+ = \{Vl + Vr), with VliVr being the spectra of the left 
and right polarisations of the transverse components and Vl = Tr since the model is parity conserving. Because of 
the condition rh < -ff* , the particle production process is found to be more efficient in the longitudinal direction than 
the transverse ones. Hence, in this case, the contribution of the vector field to the curvature perturbation is strongly 
anisotropic. Indeed, the anisotropy parameter in Eq. (|101j) in the spectrum is found to be 

9 = z(—)\ (240) 
\ m J 

where £ is defined in Eq. (|106[) and quantifies the level of the contribution of the vector field to the overall curvature 
perturbation £ of the Universe. Ref. [l2[ suggests that g < 0.3, which implies that £ <C 1 in this case. Hence, the 
vector field contribution to £ has to be subdominant, with the dominant component due to some other source, e.g. 
the inflaton field. 

Non-Gaussianity in the curvature perturbation is found to be also statistically anisotropic with its magnitude and 
direction correlated with statistical anisotropy in the spectrum as in Ref. [3f| . Indeed, for the isotropic part we found 
[cf. Eq. (HMD] 

f^l^(Z-\ = S- (241) 



ZD. A V3ff*7 3fi 

which can be substantial even if £ <C 1, where Sldcc is the density parameter at the decay of the vector field. The 
anisotropic part of /nl depends on the configuration but we have shown that Q 3> 1 in both the equilateral and local 
cases, where Q is defined in Eq. (|120p . Thus, we see that the bispectrum is predominantly anisotropic, which suggests 
that, if non-Gaussianity is observed without significant angular modulation, then the case rh < ff* in our model will 
be ruled out. 
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Consider now that rh ^> H*, which is possible only for a = —4. If this is the case, then before the end of inflation 
the mode functions begin oscillating with frequency increasingly larger than , while for the average power spectra 
we have [cf. Eq. ([ST]) ] 

1 / \ / 



Since particle production in this case is isotropic there is no statistical anisotropy in the spectrum and the bispectrum, 
i.e. g = Q = 0. Therefore, the vector field can alone generate the curvature perturbation of the Universe without any 
contribution from other sources. The model can indeed generate non-Gaussianity in the curvature perturbation with 

fi*L = -£— , (243) 

where Qa is defined in Eq. (|110[) . The above result is identical to the scalar curvaton mechanism. If there arc other 
significant contributions to £ beyond the vector field then /nl is smaller than the above. 

It is interesting to consider the case when rh ~ H*, again possible only when a = —4 (i.e. / oc a~ 4 ). In this case 
the mode functions are about to begin oscillating at the end of inflation, while their value is comparable but not 
necessarily identical. This case allows the possibility that \8V\ < V+, where 5V = V\\ — V+. Then, the vector field 
can alone generate the total curvature perturbation, but it can also produce statistical anisotropy in the spectrum 
and bispectrum, which would be at equal level 

SV 

g~g = — < 0.3. (244) 
I + 

In this case we see that the anisotropy in the bispectrum is subdominant, with the dominant component given 
by Eq. 

To investigate the parameter space for this model we have studied the evolution of the zero-mode of the vector 
field assuming energy equipartition at the onset of inflation. When a = 2 we found that W ~ constant throughout 
inflation. In contrast, when a = —4 we found that the zero mode during inflation scales as W oc a~ 3 during the 
power-law regime, while its amplitude scales as ||W|| oc a~ 3 during the oscillations, when the vector field becomes 
heavy. In all cases though, the density of the vector field pA remains constant during inflation, regardless whether 
the field is oscillating or not. After inflation, the density scales as radiation pa oc a -4 or matter pa oc a -3 when the 
field is light rh < H(t) or heavy rh > H(t) respectively. This is different from scalar fields, whose density remains 
constant when they are light. We have also verified that, when the field is heavy and undergoes oscillations, it acts 
as a pressureless isotropic fluid, which can dominate the Universe without generating a large-scale anisotropy, in 
accordance to the findings in Ref. [l[ . 

First, we considered the case when rh 3> -ff*, possible only if a = —4. As there is no anisotropy in this case we 
assumed that the vector curvaton alone generates £. By taking into account all relevant bounds on the decay rate 
of the inflaton and vector curvaton fields we have found that the scenario works when if* > 10 9 GeV. The vector 
curvaton begins its oscillations before the end of inflation but no earlier than in the last few e-folds, as N osc < 4. 
Still, this allows an exponentially large parameter space for the value of rh, which may be as large as 10 6 -H*. The 
parameter space is reduced if the decay of the inflaton is late or if the contribution of the vector field to the energy 
budget is very small during inflation. This, however, can be counteracted if the vector field decay rate is also small. If 
rh ~ then the vector field can alone generate C but may also produce statistical anisotropy within the observational 
bounds. In this case it is easy to show that -ff* > 10 9 GeV as well. Finally, the case when rh < H* allows a slightly 
lower inflationary scale since the lower bound to is relaxed by a factor < 10 _2 v / g. In this case the vector curvaton 
contribution to the curvature perturbation is strongly anisotropic, which means that Qa < 1 at decay. 



XI. CONCLUSIONS 



In this paper we studied a particularly promising vector curvaton model consisting of a massive Abelian vector 
field, with a Maxwell type kinetic term and with varying kinetic function / and mass m during inflation. The model 
is rather generic, it does not suffer from instabilities such as ghosts and may be naturally realised in the context of 
theories beyond the standard model such as supergravity and superstrings. 

We have parametrised the time dependence of the kinetic function as / oc a a , where a = a(t) is the scale factor. 
Our model offers two distinct possibilities. If rh < (possible for a = — 1 ± 3) the vector field can only produce 
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a subdominant contribution to the curvature perturbation £, but it can be the source of statistical anisotropy in 
the spectrum and bispectrum. In fact, non-Gaussianity in this case is predominantly anisotropic, which means that, 
if a non-zero /nl is observed without angular modulation, then our model is falsified in the m < case. The 
second possibility (possible for a = —4 only) corresponds to rh > H*. In this case the vector field can alone generate 
the curvature perturbation £ without any contribution from other sources such as scalar fields. If m> iJ* particle 
production is isotropic and the model does not generate any statistical anisotropy. The vector field begins oscillating 
a few e-folds before the end of inflation but its density remains constant until inflation ends. The parameter space for 
this case can be exponentially large, i.e. 1 <C m/H* < 10 6 . Significant non-Gaussianity can be generated, provided 
the vector field decays before it dominates the Universe, in which case /nl is found to be identical to the scalar 
curvaton scenario. In other words, if m> H*, our vector curvaton can reproduce the results of the scalar curvaton 
paradigm. Finally, if rh ~ -ff* the vector field can alone generate the curvature perturbation C, but it can also generate 
statistical anisotropy in the spectrum and bispectrum. In this case, the anisotropy in /nl is subdominant and equal 
to the statistical anisotropy in the spectrum, which is a characteristic signature of this possibility. 

We have also found that inflation has to occur at relatively high energies, with i?* > 10 9 GeV in the (almost) 
isotropic and > ^/glO 7 GeV in the anisotropic case (with rh > 10 TeV). These bounds correspond to prompt 

reheating, with T ro h ~ v}^ which could result in gravitino overproduction. However, if the vector curvaton dominates 
the Universe, its decay could release enough entropy to efficiently dilute the density of the gravitinos. Furthermore, 
as m Rcf. nj, one could substantially reduce the inflationary scale through introducing an increment to the mass of 
the vector field, say at a phase transition, after the end of inflation, following the mechanism first suggested for the 
scalar curvaton scenario in Rcf. [45| . 

For our model to work / and m should vary in a specific manner, which requires tuning. We have outlined two 
possibilities for achieving the desired modulation for these quantities. First, we considered that the quantities in 
question are modulated by a string modulus field, which could also play the role of the inflaton. In this case we found 
that the potential of the modulus has to be approximately exponential, which is reasonable. The second possibility 
which we discussed was that of a Higgsed vector curvaton in the context of supcrgravity theories, where scalar fields 
during inflation obtain an effective mass of order the Hubble scale. We showed that a Higgs field with tachyonic mass 
2.H* suffices to account for the desired modulation for m, while it is natural to expect that the gauge kinetic function is 
modulated by the fast-rolling scalar fields of the theory such that /// ~ H* as required. These examples demonstrate 
that the tuning of the modulation of / and m can be attained in a realistic manner. This should be contrasted with 
the traditional case of generating the curvature perturbation using scalar fields, where their effective mass needs to 
be fine-tuned at least by 0(1O -2 ) against Kahlcr corrections to produce an approximately scale-invariant spectrum 
(the famous ^-problem). In our model, the vector field also needs to be effectively massless (Af* < H*) when the 
cosmological scales exit the horizon but, to our knowledge, there is no compelling reason why this should not be 
so. Note also, that the vector field can become heavy by the end of inflation in the a = —4 case. We considered 
de Sitter inflation in our treatment so that the obtained spectra were exactly scale invariant if / cx a~ 1±3 and m oc a. 
Deviations from the Harrison- Zel'dovich spectrum can be attained either by perturbing the modulation of / and m 
or by considering quasi-de Sitter inflation, with e = —H/H 2 ~ 1CP 2 . 

Even though in our specific examples we employed scalar fields to modulate / and m, in principle their variation can 
be controlled by any kind of degree of freedom which varies during inflation. In that sense, this model can generate the 
curvature perturbation in the Universe without direct involvement of fundamental scalar fields. Of course, theories 
beyond the standard model are abundant with scalar fields and vector fields alike. The next step, therefore, is to 
realise this vector curvaton model in the context of realistic extensions of the standard model. 
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